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ABSTRACT 

A seismo-acoustic sonar concept that uses guided interface waves (Rayleigh or 
Scholte) is being developed to detect buried ordnance in the sea floor and beach 
sediments. This thesis describes the initial research conducted into the design, 
construction, and field testing of possible seismic sources that excite preferentially the 
interface waves desired for use in such a system. The theory of elasticity shows that 
seismic interface waves have elliptical particle velocity orbits in the vertical plane along 
the path of propagation. It was therefore decided that to selectively excite the desired 
interface waves, a harmonic source employed at the interface must induce elliptical 
particle motion in this plane. Several exploratory sources were developed to produce this 
type of excitation. Field tests of the discrete-mode sources developed were conducted to 
evaluate this hypothesis, but due to the non-optimum nature of the experimental sources, 
perfect discrete source excitation was not obtained. However, it was found that the 
medium itself acted as a selective filter for the interface waves after a few tens of meters 
of propagation. The experimental results obtained here suggest that the basic concept of 
discrete-mode excitation looks promising. 
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I. 



INTRODUCTION 



The research described in this thesis is a continuation of work started at the 
Applied Research Laboratories of the University of Texas at Austin, sponsored by the 
Office of Naval Research, into the development of a seismo-acoustic sonar for the 
detection of buried ordnance using guided, seismic interface waves. Their research 
identified certain key issues requiring resolution prior to further concept development. 
The primary focus of the present research is to address and resolve some these issues, 
specifically; the design, construction, and field testing of a discrete-mode excitation 
source for seismic interface waves. 

This thesis contains six chapters dedicated to this purpose. The introduction 
describes the necessity and motivation for this research. It also includes a brief 
exposition on current and emerging mine warfare technology as it relates to buried 
ordnance detection. The second chapter provides an introduction into general seismology 
with a detailed description of interface waves and the characteristics that make them 
ideally suited for such an application. The third chapter describes the previous research 
efforts conducted in this and similar areas. The following two chapters contain the major 
portion of the information on the discrete-mode source development and the experimental 
results, through the ground and surf zone testing phases. The final chapter provides 
concluding remarks and recommendations for continuing research efforts, which should 
serve to further future project development. 

A. RESEARCH OBJECTIVE 

In November of 1996, during the Mine Warfare Symposium held at the Naval 
Postgraduate School in Monterey, California, the concept of a seismo-acoustic sonar for 
the detection of buried ordnance was introduced to the mine warfare community, and was 
based on the work conducted by the Applied Research Laboratories of the University of 
Texas at Austin (ARL:UT). [Ref. 1] 

The main objective of the ARL:UT study was to demonstrate that an electro- 
mechanical transducer could be employed to excite seismic surface waves, of which there 



are many types, but especially Rayleigh waves, which could be used to ensonify a buried 
mine-like target, and which in turn would produce seismic echoes that could be received, 
processed, and displayed. Due to the generally large contrast in density and elasticity 
between buried ordnance and the surrounding soil, the target produces seismo-acoustic 
echoes that travel in all directions, in all of the possible propagating modes allowed by 
the theory of elasticity [Ref. 1]. The experimental results of the study showed that the 
concept of a seismo-acoustic sonar for the detection of buried ordnance was scientifically 
feasible. 

During the course of this first investigation into the use of a seismo-acoustic sonar 
for buried ordnance detection, certain key issues requiring resolution were identified. 
These issues need to be solved before the overall concept will have sufficient credibility 
to progress onto naval system development. These issues included [Ref. 1]: 

• Development of sources for selective, or discrete-mode, excitation of seismic 
interface waves. 

• Development of techniques for selective reception of seismic target echoes. 

• Adaptive methods for excitation and reception of optimal modes. 

• Measurement, analysis, and utilization of real mine target strengths. 

• Incorporation of new discovered cetacean sonar protocols and signal physics. 

• Re-examination of vector signal processing technology. 

The purpose of the present research is to address and attempt to resolve some of 
these key issues, specifically, to initiate development of a discrete-mode source and to 
further refine the vector signal processing techniques. The results of this research should 
bring the concept of a seismo-acoustic sonar closer to field implementation. 

B. RESEARCH MOTIVATION 

Because of the series of events that transpired during the Persian Gulf War, a 
great deal of attention has been directed towards “main-streaming” mine warfare in the 
US military. The use of mines in the battlefield poses a serious threat to many different 
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military operations, and thus significantly impacts the United States’ ability to execute its 
national security objectives. The tactical use of mine warfare directly influences the 
shape of the battlespace and may exclude our military forces from areas in which they 
our called upon to serve. In order to accomplish its military objectives, US forces need 
effective methods to detect and counter enemy mines. As Admiral Boorda stated in his 
1995, white paper on mine warfare, “This warfare area must take its proper place among, 
as well as influence, future research and development, acquisition, and modernization 
plans and programs’’[Ref. 2]. 

In terms of cost-effectiveness and operational and psychological impact, mines 
are perhaps the single most attractive weapons available for disrupting the battlespace. 
For the foreseeable future, we must expect increases in both the lethality and number of 
mines available. Therefore, more advanced technologies must be pursued to make the 
detection, classification, and removal of mines and other buried ordnance items quicker, 
safer, and more cost-effective. 

C. CURRENT AND EMERGING MINE WARFARE TECHNOLOGY 

At the present time (1998) there is no single method currently available that can 
detect buried mines in shallow water sediments, through the surf zone, onto the beach and 
beyond. All the varying buried mine amphibious countermeasure operations can be 
effectively grouped into either breaching for amphibious assaults or post conflict 
clearance operations. Breaching operations involve the construction of safe paths 
through a minefield during combat, while clearance operations deal with eliminating 
unexploded ordnance from large areas for occupation and use. Accurate detection and 
identification are essential for both types of operations. For combat breaching it is 
important to quickly know where the minefields are located and what their density is. 
while ordnance clearance operations need to deal with each and every mine, on a non- 
real-time basis. 
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1 . 



Current Mine Warfare Technology 



a) Shallow Water Mine Warfare Technology 

From the high-water mark seaward Explosive Ordnance Disposal (EOD) 
divers, Naval Special Warfare (NSW) combat swimmers, Marine Corps reconnaissance 
swimmers, and the U.S. Navy marine mammals, plus a variety of towed sweep gear 
provide the bulk of the available shallow water mine countermeasure systems. 

EOD and NSW personnel serve as an extremely effective neutralization 
component, but provide only minimal detection capability. Sonar and laser detection 
systems such as the Coastal Battlefield Reconnaissance and Analysis system (COBRA), 
the Remote Minehunting System (RMS), and the Airborne Laser Mine Detection System 
(ALMDS) are limited to detection of floating mines, shallow water moored mines, and 
exposed bottom mines. Of all the sonar and laser technologies currently employed for 
shallow water detection only one system is truly effective in the detection of buried 
mines, namely the marine mammals. Even though great strides have been made in 
increasing the effectiveness of marine mammal systems, search rates are low and the 
mammals have difficulty working in the surf zone or in warm water. 

b) Land Warfare Technology 

From the high-water mark landward, the currently employed breaching 
and clearance technologies revolve around hand-held probes, metal detectors, ground 
penetrating radar, trained dogs, and mechanical breaching equipment. Each of these 
methods has its unique advantages and disadvantages, but they are in general slow, 
costly, labor intensive, and not optimal for combat operations. 

Manual mine countermeasures provide the most reliable and consistent 
method of ordnance location, classification, and identification. But manual mine 
clearance is slow and dangerous, with an injury rate of nearly one technician every 1,000 
mines cleared [Ref. 3j. 
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“Metal detectors” are only capable of detecting ordnance containing 
ferrous components. They are generally difficult to successfully operate because they 
have high false alarm rates. The accuracy of these systems to detect, locate, and analyze 
target returns is based on operator judgment and experience. The high false alarm rates 
associated with these systems is usually due to the high ferromagnetic soil content in 
battlefield areas, as well as by the presence of man made ferromagnetic false targets such 
as hubcaps, tea pots. etc. [Ref. 4] 

Ground penetrating radar incorporates the same general technology as 
conventional radar, but on a much smaller scale. Ground penetrating radar systems fall 
into four main categories: common short pulse systems, video pulse systems, synthetic 
aperture pulse radar, and frequency modulated-continuous wave radar [Ref. 4], They all 
use the propagation of electromagnetic energy into the soil for target detection and 
localization. Subsurface sensing radars rely on the strength of the electromagnetic signal 
reflected from the target. The strength of the reflected signal is highly dependent on 
target distance, composition, and orientation, in addition to soil characteristics such as 
layering, presence of boulders, etc. High resolution ground penetrating radar systems 
usually operate in the frequency range of 30 MHz to 2 GHz, which implies short 
operating ranges due to the high attenuation of electromagnetic energy in the ground at 
these frequencies. They also suffer from high false alarm rates. 

Trained dogs provide a high degree of detection capability without the 
requirement for prior visual location. Most explosive sniffing dogs are trained to identify 
only one or two types of explosives; thus a whole squad of dogs would be needed to 
detect all the possible explosives types used in just landmines alone. Additionally, the 
amount of time required to train a dog for buried ordnance detection is extremely high, 
while search rates are low because the dogs get tired and detection capability is reliable 
only over very short ranges. 

Lastly, breaching equipment is the brute force method of land mine 
countermeasures. The use of breaching equipment is the most expedient, but of the low- 
end types of available breaching equipment, rollers rarely achieve more that 70% 
clearance effectiveness and rotary flails or blades less than 80%. Explosive net breaching 
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whereby a “fish net'’ containing shaped charges is rocket fired over an area is more 
effective, but much more expensive. 

Although these current methods offer the greatest assurance that an area is 
safe for use, they are also quite dangerous because they put the operator in close 
proximity to the ordnance. 

2. Emerging Mine Warfare Technology 

Emerging technologies with application to detection and clearance can be grouped 
into the following categories: infrared, microwave, photon (laser linescan) backscatter, 
nuclear or thermal neutron analysis, passive optical imaging, and “new” acoustic and 
seismic technologies. Many of these approaches suffer from large power and high data 
processing requirements, and their effectiveness is highly dependent on the propagation 
characteristics of the medium. 

D. PROPAGATION MEDIUM CHARACTERISTICS 

To realistically address the issue of buried mine detection in the surf zone and on 
land, much emphasis must be placed on understanding the characteristics of the media in 
which the mine detection system must operate. Depth of signal penetration, resolution 
characteristics, detection capability, and other system parameters are all highly dependent 
upon the medium of propagation. These conditions will usually vary from location to 
location and from season to season. 

1. Soil Characteristics 

In general, the earth’s boundary layers consist of a variety of minerals, rocks, 
organic and inorganic compounds, water, and air. The amount of organic materials 
present is specific to particular areas and geographical regions, but most organic 
materials significant to mine warfare are only present in the first few meters of the upper 
sediment layer, and decrease with increasing depth. Changes in the water content result 
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in changes of the dielectric and conduction characteristics of the medium. All these 
varying conditions have significant effects on the propagation of electromagnetic and 
acoustic waves. [Ref. 4] 

Nominally, low frequency electromagnetic waves are less attenuated and result in 
increased depths of penetration. Higher frequency electromagnetic signals are attenuated 
more, but provide better resolution characteristics. Acoustic energy is less attenuated 
than electromagnetic energy in high moisture and densely compacted soil environments. 

2. Surf Zone Characteristics 

In offshore sediments and in the surf zone, the affect of the water and sediment 
interfaces on energy propagation are quite pronounced. Optical techniques have been 
developed to penetrate the water interface, but cannot penetrate the underlying sediment, 
so are essentially useless in buried ordnance detection. Magnetic techniques are also 
limited because of their low resolution qualities, short operating ranges, and the increased 
use of non-magnetic mines. Radar systems are ineffective in the surf zone because of the 
extremely high attenuation suffered by electromagnetic energy in seawater, penetrating 
no more than the first meter into the water. Acoustic sonar systems are capable of 
penetrating both the water and sediment boundaries, but in the surf zone extensive clouds 
of micro bubbles generated by wave action scatter and attenuate most of the acoustic 
energy. The energy which is left can penetrate the sediments, but is limited by Snell's 
Law, which establishes a “critical grazing angle’’, below which sound is totally reflected 
within the water column and does not penetrate the sediments. [Ref. 5] 

E. SEISMIC WAVE RESEARCH 

With the limitations imposed by the operating environment to current buried 
ordnance detection systems, it was obvious that another method of energy propagation 
must be identified that could more effectively utilize the soil and sediments as a 
propagating medium. With this in mind, it became necessary to extend research efforts to 
incorporate seismic wave techniques. Seismic surface waves also propagate along both 
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water/sediment and air/sediment interfaces and can be easily excited in the environments 
of primary interest for buried mine detection systems. Because of their propagating 
characteristics, the use of seismic interface waves can provide a safe minimum standoff 
distance between the operator and the ordnance. 

There are no systems currently available today that can detect buried mines on 
land and in the surf zone. It is clear that development of such a system would help to 
provide the military with the kind of system needed to make the transitions to the shore, 
imp. ;d in such documents like “...From the Sea”, physically realizable. The benefits to 
other operations, such as, humanitarian demining could also be significant. 

The most common types of seismic waves are those generated by earthquakes, 
which travel both through the body of the earth and along its surface. Seismic waves 
have also been used extensively in geological prospecting for oil and mineral deposits. 
The next chapter contains an overview of seismic wave theory and specific information 
on interface waves and their application to buried ordnance detection. 
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II. SEISMIC WAVE THEORY 



This chapter contains a description of the types of seismic waves that propagate 
through an elastic medium and along its surface. Within this cornucopia of wave types 
exist a category of boundary waves known as Rayleigh and Scholte waves. The general 
characteristics exhibited by these waves that make them of practical use in buried 
ordnance detection systems will be explored as well as their relation to the other seismic 
waves. 



A. GENERAL SEISMOLOGY 

When one studies the general theory of elasticity, two types of mechanical waves 
can be found to propagate in an unbounded, infinite, homogeneous, elastic medium. 
These waves, one compressional and one shear, are called body waves since they travel 
within the body of the propagating medium. In 1885, the English scientist. Lord 
Rayleigh (J.W. Strutt) [Ref. 6] demonstrated theoretically that there exist waves that can 
propagate over the planar interface between a semi-infinite elastic half-space and a 
rarefied gas. These waves later came to be known as Rayleigh waves. Further research 
lead to the discovery of many other boundary waves. Seismic waves can be grouped into 
one of the following two categories: body waves and surface waves. 



1. Body Waves 



There are two types of body waves. The first and fastest propagating wave is 
called the Primary wave or P-wave. A P-wave is a compressional wave, in which the 
wave motion alternately compresses and dilates the medium. The particles of the 
medium undergo displacements in a direction co-linear with the direction of propagation: 
hence the P-wave is also called the longitudinal wave. P-waves can travel through both 
solid rock and liquid materials with a (longitudinal) propagation speed given by 



P 






( 2 . 1 ) 
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where B is the adiabatic bulk modulus, p the shear modulus, and p the density of the 
medium. 

The second, and slower, body is wave is called the Secondary wave or S-wave. S- 
waves are shear waves, in that, as an S-wave propagates it shears the medium at right 
angles to the direction of propagation. Due to the fact that liquids do not support finite 
shear stresses, S- waves cannot propagate in liquids. S- waves travel with a speed given 
by 




In general, when an external force deforms a solid, both P-waves and S-waves 
will normally be produced. When a wave of either type interacts with an interface 
boundary, reflected waves of both types can be produced. 

2. Surface Waves 

When a solid has a fluid surface (air or water), or when there is an interface 
boundary between two media or between layers of a single medium, surface waves are 
normally present. Surface waves are somewhat similar to gravitational ocean waves for 
which most of the wave motion is located near the surface and decays exponentially with 
increasing distance into both the elastic substrate and the overlying fluid, although the 
mechanisms are different. 

Surface waves can also be divided into many types. One type of surface wave, 
which propagates in the duct between two sedimentary boundaries, is called a Love 
wave. Love waves induce side-to-side particle motion in a horizontal plane, but at right 
angles to the direction of propagation. A second type and more prevalent surface wave is 
the Rayleigh wave. Rayleigh waves, which occur at the planar interface between a semi- 
infinite elastic half-space and a gas, have both vertical and horizontal particle motion in a 
vertical plane oriented along the direction of propagation. Rayleigh waves propagate 
with a phase velocity slightly less than the bulk shear wave velocity in the elastic 
medium. Another type of boundary wave, analogous to the Rayleigh wave, is the Scholte 
wave, which occurs at the interface between a liquid and an elastic half-space such as the 
seafloor. Scholte waves are characterized by a phase velocity that is also slightly slower 
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than the shear wave speed in the solid. Essentially, Rayleigh and Scholte waves are very 
similar, the only difference being that the liquid media component of Scholte waves is 
more pronounced than the gas media component of Rayleigh waves. The interface wave 
that propagates along the boundary created by two elastic half-spaces is called a Stoneley 
wave. Stoneley waves travel with a velocity that is slower than the shear wave velocity 
in either media. 



B. RAYLEIGH WAVES 



In an isotropic solid the stress-strain relations have the following form 
+ 2pe„ , cr v> , = >-0 + 2\x.z yy , a„ = >.0 + 2p e :: ^ 

a xy = , ct v , = pe*. , a., = pe„ 

where X and p are the Lame constants, 0 the dilatation, or volume strain, and <j y and e y 
are the stress and strain, respectively, where / denotes the direction of the stress or strain 
and / defines the plane over which it is applied. Using the elastic relations from (2.3), the 
equations of motion for an isotropic solid (neglecting body forces) can be written as 

p0=(A. + p)ra + nF*«, (2.4) 



where u represents the particle displacement and p the density of the medium [Ref. 7]. 
These equations of motion can be shown to lead to the propagation of the two types of 
body waves with phase velocities given by (2.1) and (2.2). 

If we consider the propagation of a plane wave through an elastic medium with a 
free boundary, such that the boundary is free from stress, we are lead to the wave 
equations 



d t 2 



= cf ?< )>, 



d Vg 

dt 2 



c; lAy , 



(2.5) 



where the potentials <j> and \\i were introduced to separate the effects of compression and 
rotation in the medium [Ref. 7]. The simultaneous solution of (2.5) for a disturbance that 
is localized at the boundary and decays with increasing depth requires that the following 
relation be satisfied 

k* - SiU + (24-1 &i*y* +(l6y 2 -16) = 0. (2.6) 
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This relation is known as the Rayleigh equation, where k is the ratio of the Rayleigh 
wave speed to the shear wave speed, normally 0.9, and y 2 is the ratio of the shear wave 
speed to the longitudinal wave speed, y 2 can also be expressed solely in terms of the 
Poisson’s ratio, namely 



t=t = 



p 



1 -2v 



'/ X + 2p 2 - 2v ' 



(2.7) 



where v is the Poisson’s ratio. Thus, to solve the Rayleigh equation, which is cubic in k 2 , 
requires that only the Poisson’s ratio for the propagating medium be known. Therefore 
the Rayleigh wave speed in an infinite isotropic solid is independent of frequency and 
depends only on the elastic constants of the medium. The Rayleigh wave corresponds to 
the root of the Rayleigh equation that lies in the (real) interval (0,1). It can be shown that 
for values of the Poisson’s ratio in the range 0 - 0.5 the Rayleigh equation has one and 
only one such root. Hence Rayleigh waves (R-waves) traveling at the boundary of an 
isotropic, homogeneous medium are dispersionless. [Ref. 7,8] 

Figure 2.1 shows the variation in the wave speed ratios as a function of the 
Poisson’s ratio for an isotropic, homogeneous solid substrate. From the figure we can see 
that Rayleigh waves propagate at almost the same speed as shear waves, varying from 




Figure 2.1 . Variation in Wave Speed Ratios as a Function of Poisson’s Ratio 
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approximately 0.87 to 0.96 of the shear wave speed. Appendix A1 contains the 
MATLAB code used to generate Figure 2.1 and Appendix A2, the code used to evaluate 
just the Rayleigh wave speed ratio as a function of the Poisson’s ratio. 

Rayleigh waves can be visualized as the superposition of two waves, one 
longitudinal and one transverse, which propagate along the boundary of an elastic half- 
space. The particle motion associated with a Rayleigh wave is normally elliptical in 
nature, with the semi-major axis of the ellipse perpendicular to the boundary and the 
semi-minor axis parallel to the direction of propagation, depending on the boundary 
conditions. As Figure 2.2 illustrates the particle motion at the surface is retrograde 
elliptical, but shifts to prograde at a finite depth. The shift in particle rotation occurs at a 
depth of approximately 0.2X r for solids, and approximately 0.1 X r for real sediments, 
where X r is the Rayleigh wavelength [Ref. 8,9], 




Figure 2.2. Elliptical Particle Motion Induced by Rayleigh Wave 

Figure 2.3 displays the variation in the ratio of the longitudinal to vertical particle 
displacements for Rayleigh waves at the surface of an homogeneous, isotropic solid, as a 
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normalized displacement amplitude longtudnal dsdaoemertAetical dsplacemsr* 



function of the Poisson’s ratio, while Figure 2.4 depicts how the particle displacements 
vary with depth. 




Figure 2.3. Ratio of Longitudinal to Vertical Particle Displacements 
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Figure 2.4. Variation of Longitudinal and Vertical Particle Displacements 
(normalized by the surface vertical particle displacement ) 
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Figure 2.5 shows the variation in the stress components that occur with depth for the 
same isotropic conditions. 




Figure 2.5. Tensile and Shear Stresses 



From Figure 2.4 we can see that the vertical displacement first increases and 
reaches its maximum at approximately 0. 1 ~k r , then exponentially decreases with depth. 
The longitudinal displacement decreases with depth from the outset and changes 
direction at approximately 02X r , reaching a secondary maximum at 0.5A. r . From Figure 
2.5 one can see that the tensile stress along the direction of propagation also changes 
direction at approximately 025X r and the other stress components reach their maximums 
at 0.3 X r , then all exponentially decay with depth. These sets of curves serve to show the 
localization of the Rayleigh wave to a layer approximately two wavelengths thick under 
isotropic conditions. 

Appendices A3-A6 contains the MATLAB code used to generate Figures 2.2-2. 5, 
respectively. Additional insight into Rayleigh wave behavior is contained in the 
computer codes and figures making up the rest of Appendix A. 
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c. 



SEISMIC WAVES 



From the literature [Ref. 6,7,10], where solutions to the elastic wave equations 
have been obtained, we can form an idealized picture, as shown in Figure 2.6, of the 
wave system created by a vertical point impulse at the surface of an ideal, isotropic. 





Figure 2.6. Wave System from Vertical Point Impulse [Ref. 10] 

substrate medium. From the figure we can see that the wave system is composed of 
three major events corresponding chronologically to the arrivals of the P-wave, bulk S- 
wave, and R-wave. Thus a particle at the surface will first be displaced by the P-wave, 
then by the bulk S-wave, followed by a much larger oscillation due to the passing 
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Rayleigh wave. This illustrates that the Rayleigh wave is the most significant disturbance 
along the surface of an isotropic half- space, and at large distances from the source may 
be the only clearly distinguishable wave. [Ref. 10] 

For a disturbance in the form of a circular footing undergoing vertical oscillations 
at the surface of an isotropic, elastic half-space the basic features of the wave field far 
away from the source are shown in Figure 2.7 [Ref. 10]. 



Circular Footing 
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Figure 2.7. Wave Field from a Circular Footing 



The two body waves propagate radially outward from the source along a hemispherical 
wavefront. The decrease in the amplitudes of the body waves within the medium is 
proportional to 1/r, where r is the radial distance from the source, while along the surface 
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of the half-space the far-field amplitudes decrease as 1/r 2 . Rayleigh waves, on the other 
hand, propagate radially outward with a cylindrical wavefront and experience a far-field 
amplitude decrease that is proportional to 1/Jr. For a vertically oscillating point source 
at the surface of an isotropic, elastic half-space, it was determined that the distribution of 
input energy among the three wave types is 67% Rayleigh wave, 26% shear wave, and 
7% compressional wave. [Ref. 10] 

From this brief exposition on seismic wave fundamentals, the features that make 
Rayleigh waves unique (in an isotropic half-space) are that they are: 

• Dispersionless, thus propagate with a phase velocity that is independent of 
frequency. 

• Localized within a layer that is approximately two wavelengths deep. 

• The most significant disturbance along a surface boundary. 

• Easy to excite, with a particle motion that is retrograde elliptical at the 
surface. 

It is these features that make Rayleigh waves ideal for buried ordnance detection 
in the earth and along the beach. Scholte waves exhibit the same general characteristics 
as Rayleigh waves, the major difference being that Scholte waves occur along the 
interface between a liquid and an elastic half-space, thus making them ideally suited for 
surf zone applications. 

It should be remarked that in natural soil, beach, or seafloor sediments, the 
overburden weight of the surface sediments causes the shear modulus to increase with 
increasing depth. This means that the shear velocity increases with depth and that, in 
turn, means we have a condition requiring dispersion. This occurs because the lower 
frequency, long wavelength waves naturally propagate at deeper depths in the sediment, 
where the shear velocity is greater. Thus, low frequency, long wavelength signals 
propagate faster than high frequency, short wavelength signals. 

If an electromechanical source could be built that could selectively excite 
Rayleigh waves, more of the energy used in signal generation could be put on target, and 
less of the projected energy would exist to create unwanted reverberation carried in non- 
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Rayleigh wave modes. However, at the target, even pure Rayleigh waves get distributed 
into other non-Rayleigh modes. With the addition of special signal processing techniques 
that exploit the known features of Rayleigh waves from target echoes, we could capitalize 
on discrete-mode sonar methods. The next chapter discusses the previous research 
efforts made in buried ordnance detection using seismic waves and the experimental 
limitations imposed by the systems used. 
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III. 



PREVIOUS RESEARCH EFFORTS 



This chapter briefly reviews two previous research projects that explored the 
propagation of energy along an air/sediment interface to detect buried objects. The first 
project was an experiment conducted at the Naval Postgraduate School, investigating the 
applicability of array theory to target localization. The second project, done at the 
Applied Research Laboratories of the University of Texas at Austin (ARL:UT), for which 
this thesis is a direct continuation, developed the theory and signal processing techniques 
for the first phase of the seismo-acoustic sonar system. Other significant prior research 
consists of the dry-land experimental work done at the firm of Bolt, Beranek, and 
Newman (BBN) and theoretical work done at the Coastal Systems Station. These works 
will be cited subsequently. 

A. NAVAL POSTGRADUATE SCHOOL (NPS) 

In 1995, NPS student LT William Stewart conducted research into the feasibility 
of using a phased array to generate surface waves to detect buried ordnance [Ref. 1 1], His 
research was the continuation of a study conducted by BBN Systems and Technology 
from 1987 to 1992 [Ref. 12]. The source utilized in the NPS study was a four-inch, large 
cone excursion, electrodynamic speaker, whose frequency response for the radiation of 
sound was relatively constant over the range 50 - 7000 Hz. The speaker cone was 
modified by bonding to it a 25-gram mass, cone-tipped plunger, which was used to 
improve ground coupling. See Figure 3.1. 

Three single-axis seismometers were used as ground vibration receivers. Testing 
was conducted in a twelve-foot diameter, four- foot deep redwood tank filled with 
approximately 17 tons of medium grain sand. With a single source and the three 
receivers, LT Stewart was able to successfully demonstrate target localization using 
surface waves. Though no data was reported on the results for target localization using a 
phased array, his study did show that beamforming with a three element linear array of 
surface sources was possible. [Ref. 1 1] 
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B. ARL:UT 

An investigation into the use of Rayleigh waves for buried ordnance detection 
was conducted by ARL:UT on the beach of Mustang Island, near Corpus Christi, Texas, 
in 1996 [Ref. 1], Their experimental setup consisted of a source which employed a six- 
inch by eight-inch exciter foot, from which protruded 43, 3 /4-inch nails and a receiving 
array composed of three seismometers, as shown in Figure 3.2. 

Their raw, field recorded data did not show the target. But with the use some 
sophisticated vector signal processing techniques and some coherent subtraction methods, 
operating on the reverberation, the target was detected at its appropriate range and 
bearing. A seismo-acoustic detection model for elastic targets buried in elastic sediments 
was also developed to support the measurements [Ref. 13]. 
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Figure 3.2. ARL:UT Seismic Sonar Experiment 



C. COASTAL SYSTEMS STATION 

Some of the benchmark research efforts recently produced by the laboratories of 
the Coastal Systems Station concern the theoretical solutions and numerical modeling of 
acoustic wave scattering from elastic objects buried in underwater sediments and plane- 
stratified poroelastic media [Ref. 14,15]. Objects located in fluid mediums produce 
resonant scattering behaviors. These same objects when buried in ocean sediments 
display a resonant scattering behavior which is modified by the proximity to interfaces, 
sediment attenuation, and sediment loading, thus information obtained from the scattering 
field can be used for identification as well as detection. Further experimental verification 
of these solutions could have significant impact on the capabilities of seismo-acoustic 
sonar systems. 
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D. 



EXPERIMENTAL LIMITATIONS 



What is common to the aforementioned experiments is the use of vertical source 
excitation, which generates not only Rayleigh waves, but also undesired P- and S-waves. 
Thus, the Rayleigh wave signals of interest were contaminated by body wave signals. To 
satisfactorily localize the target, in the ARL:UT experiments, this unwanted reverberation 
had to be removed during the signal processing by coherent subtraction of the reverberant 
background field in the absence of the target from the field with the target present. 
Vector signal processing was also used for the first time, in the ARL:UT experiments, to 
bring the target echoes above the reverberation. This technique capitalizes on the known, 
90-degree phase difference between vertical and horizontal components in propagating 
Rayleigh waves. But since coherent subtraction is not a feasible signal processing option 
in a operationally fielded system, another solution must be identified. The research team 
at ARL:UT theorized that considerable reverberation suppression could be achieved by 
utilizing an excitation source that could selectively excite Rayleigh waves [Ref. 1]. They 
also experimented with a two-degree of freedom source, which preferentially excited 
Rayleigh waves that propagated to tens of meters in range. Unfortunately, this source 
produced a signal too long in time duration for sonar application. 

The use of ultrasonic Rayleigh waves in non-destructive testing is a well- 
documented science, but selective excitation of low frequency Rayleigh waves in 
terrestrial media has not been previously attempted. The remainder of this thesis 
describes the development and testing of a discrete-mode, electromechanical, Rayleigh 
wave source for a seismo-acoustic sonar. 
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IV. DISCRETE-MODE SOURCE DEVELOPMENT AND TESTING 

(GROUND PHASE) 



This chapter documents the initial design, construction, and testing of the first 
four generations of the discrete-mode, Rayleigh wave source. The associated 
experimentation involved testing in both a laboratory environment and various field 
environments. This portion of the research, known as the Ground Phase, was used for 
proof-of-concept testing and qualitative determination of the relations between input 
signal parameters and the resulting surface wave field. 

A. DESIGN CRITERIA 

The first step in constructing a discrete-mode, electromechanical, Rayleigh wave 
source for a seismo-acoustic sonar system was to establish a set of guiding design 
criteria. Based on the known features of Rayleigh waves and the previous research work 
completed, it was decided that minimum requirements that a source needed to meet were 
that it: 

• Have sufficient degrees of freedom (2) to drive the ground in elliptical orbits. 

• Have good ground coupling in both the horizontal and vertical directions. 

• Generate Rayleigh waves with sufficient amplitudes to allow detection at a 
range of at least 10 m. 

• Provide a repeatable waveform capable of high repetition rates over long 
periods of time. 

• Be able to operate over a range of low frequencies in the 20 - 200 Hz band. 

The most important of these requirements is that the transducer oscillate the 
ground in elliptical orbits, so that as much as possible of the input energy be converted to 
discretely excite Rayleigh waves. 

All seismic transducers consist of two basic components: a vibration source and a 
source-to-ground coupling mechanism. The vibration source transfers the input electrical 
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energy into output mechanical energy. The source-to-ground coupling mechanism 
provides the means by which the output mechanical energy is transmitted into the 
propagating medium. A constraint on the source design was that the ground coupling 
mechanism be as light weight as possible, while at the same time being as stiff as 
possible. This was desired so most of the mechanical force generated by the source 
would be used to move the ground and not just the source. 

B. CONCEPT SOURCE 

With the design criteria established for a suitable source, the next step was to 
design and construct laboratory prototypes that could be tested to examine the effects of 
varying source parameters such as frequency and driving voltage, on the resulting surface 
wave field. This was deemed a necessary first step because the results obtained would 
help in selecting an appropriate vibration source and designing the ground coupling 
mechanism for a field test source. It was decided to construct the concept source out of 
on-hand laboratory equipment and to conduct the required testing in the twelve-foot 
diameter, four-foot deep redwood, sand filled tank used in the previous NPS experiment. 

1. Vibration Source 

The Model F5B Electromagnetic Shaker (or Mini-Shaker for short) produced by 
Wilcoxon Research was chosen as the vibration source for proof-of-concept testing. The 
Mini-Shaker is an electrically driven force generator, with a constant (blocked) force 
output of 0.2 lbs ( 0.89 N) over the frequency range of 100 - 1000 Hz (See Appendix B1 
for manufacturer specifications). The design configuration of the Mini-Shaker allows it 
to be stud-mounted in any position without external support or shaft alignment 
considerations. 
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2 . 



Ground Coupling Mechanism 



The source-to-ground coupling mechanism, as shown in Figure 4.1, was 
constructed from two pieces of aluminum. Aluminum was chosen because of its high 
tensile strength and low mass. Appendix C 1 contains the actual design drawings. 




Figure 4. 1 . Picture of Concept Rayleigh Wave Source (scale in inches) 



The bottom piece was constructed from a hollowed out aluminum block. When 
buried in the sand, the inside surfaces of the block provide the required vertical and 
horizontal coupling to the medium. The top piece was manufactured to hold two Mini- 
Shakers and a tri-axial accelerometer (Accelerometer specifications contained in 
Appendix B2). Mounting studs were placed in the top piece so as to align the Mini- 
Shakers at a 45 degree angle with the horizontal surface plane of the sand. This specific 
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placement of the Mini-Shakers was chosen so that any two-dimensional oscillatory 
motion could be excited by individually controlling the phase and amplitude of each 
shaker. 



3. Experimental Results 

The test equipment setup as described in Appendix D1 was used for the concept 
source testing. The source and seismometer receivers were positioned in the test tank as 
shown in Figure 4.2 and Appendix Dl. This setup allowed for the greatest separation 
between the source and the tank wall and the farthest receiver and the tank wall. The 
seismometers were tri-axial, moving coil geophones, which measure the velocity of 
ground movement. 




Figure 4.2 Concept Source Testing Experimental Setup 
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The first experiment was intended to qualitatively determine the effects of phase 
differences between the Mini-Shakers on the axial surface wave field, as a function of 
frequency. These results were characterized by the difference in the seismometer outputs 
for the axial and vertical directions. The transverse direction was not evaluated because 
the Rayleigh wave motion of interest exists in the vertical plane along the direction of 
propagation. Continuous sinusoids of frequencies 100, 150, 200. 250, 300, and 500 Hz 
were chosen as the test values. Since the left most Mini-Shaker in Figure 4.2 was used as 
the reference shaker, phase differences were chosen to correspond to the following six 
cases: 

• Case 1 - Right shaker lagging by 25 degrees. 

• Case 2 - Both shakers in phase. 

• Case 3 - Right shaker leading by 25 degrees. 

• Case 4 - Right shaker leading by 45 degrees. 

• Case 5 - Right shaker leading by 90 degrees. 

• Case 6 - Right shaker leading by 1 35 degrees. 

In all test cases except the second, which produced exclusively vertical motion, the 
motion of the source was elliptical as was seen by the lissajous pattern of the radial and 
vertical components of the accelerometer. Appendix E contains the qualitative 
experimental results for all phases and frequencies evaluated. The results of this first 
experiment showed that there was a great deal of variation in the axial surface wave field 
with both source frequency and relative phase. 

A second experiment was conducted to assess the effects of not only varying the 
frequency and phase, but also the driving voltages of the Mini-Shakers. This second 
experiment further served to show that the axial surface wave field, even in the nearly 
homogeneous, isotropic conditions of the sand tank, was a highly complex outcome of all 
variable source parameters. 

Even though surface waves were generated in all test cases, the specific 
combination of source parameters needed to produce Rayleigh waves of sufficient 
amplitude and duration at the desired ranges for buried ordnance detection had yet to be 
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determined. It was evident from the initial testing that the appropriate combination 
would not be found in the laboratory. This is because, not only is the correct 
combination highly dependent on the individual properties of the vibration source and the 
ground coupling mechanism, but also on the local properties of the propagating medium 
at the time of testing. Furthermore, the limitations imposed by the tank, such as 
reflections from the tank walls and short ranges, made it impossible to quantify the 
source. To proceed any further in developing a discrete-mode, Rayleigh wave source 
would require designing a source that could be tested under actual field conditions. 

C. FIELD TEST SOURCE -MOD I 

To begin constructing a discrete-mode, Rayleigh wave, source that could be used 
in the field required that a suitable vibration source be identified. Once the vibration 
source was chosen, then the source-to-ground coupling mechanism could be 
appropriately designed and constructed. The vibration source needed to be rugged 
enough to withstand the beach and surf zone environment, yet inexpensive enough so that 
multiple pairs could be acquired, if arrays of paired vertical and horizontal sources could 
eventually be utilized. 

1. Vibration Source 

A test source that satisfied most of the necessary prerequisites was found 
manufactured by Aura Systems Inc., under the nomenclature of a Bass Shaker (See 
Appendix B3 for manufacturer specifications) [Ref. 16]. The Bass Shaker is a moving 
magnet transducer, intended to be used as a low frequency vibration source for a car 
stereo system. They are normally mounted under car seats or in door panels to provide 
listeners with the sensation of sound, without moving air like an ordinary speaker. The 
Bass Shaker generates a nominal force output of 10 lbf (44.5 N), independent of its 
mounting orientation. 
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2. Ground Coupling Mechanism 

Since the Bass Shakers were considerably different in size and shape from the 
Mini-Shakers used previously, the source-to-ground coupling mechanism had to be 
redesigned. It was decided that the Bass Shakers should still be positioned at a 45- 
degree incline to the horizontal plane to provide the necessary degrees of freedom to 
produce all variations of elliptical motion. As shown in Figure 4.3, the ground 




Figure 4.3. Picture of Field Test Source - Mod 1 



coupling mechanism was constructed from two pieces of aluminum. The top piece 
serving as a support brace for the upper part of the Bass Shakers and as a mounting 
surface for the electrical (BNC) connectors. The bottom piece serves as the attachment 
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point for the lower part of the Bass Shakers and provides the coupling mechanism to the 
medium through two rows of machine screws. Appendix C2 contains the Mod I 
technical drawings. A tri-axial accelerometer was positioned in the center of the base to 
monitor the source motion, and there was also an accelerometer mounted on the base of 
the support brace. Refer to Appendix B2 for accelerometer specifications. 

3. Experimental Results 

Before commencing any field experimentation, the new Rayleigh wave source 
was put through some operator familiarization and initial shakedown testing in the sand 
tank laboratory. Using the test equipment setup as described in Appendix D2, the same 
experiments were conducted as with the Mini-Shakers. These tests provided further 
verification that the axial surface wave field was a strong function of the input source 
parameters. 

With the knowledge gained from the laboratory, it was clear that the first goal of 
any field experimentation s! .Id be to assess the Mod I source’s ability to generate 
surface waves. To do this required finding a nearby location that was relatively free from 
any underground obstructions that would scatter the outgoing surface waves and make 
detection more difficult. The choice was made to use the NPS softball field as the first 
site for the field experiments. Progressing onto field experimentation from an equipment 
point of view proved to be quite a challenge at first, but as the research progressed, 
equipment difficulties were quickly overcome. 

Once at the softball field, the lab equipment was setup as described in Appendix 
D3. Figure 4.4 illustrates the initial experimental layout. The initial distances were 
chosen based on expected Rayleigh wave speeds on the order of 100 m/sec and the 
maximum length of electrical cable available. Due to the hardness of the ground, 
trenches for the two rows of machine screws had to be dug in order to get the base of the 
source flush with the ground. After all the equipment was setup, attempts were made to 
detect the source motion above the background noise at the two seismometers located at 
ranges of 15.2 m (50 ft) and 30.5 m (100 ft). After many unsuccessful attempts to detect 
Rayleigh waves at these ranges, the decision was made to move back to ground zero. 



32 



The two seismometers were repositioned at distances of 1 m and 5 m from the source. 
From the variation in source parameters attempted, it was clear from the accelerometer 
output that the Mod I coupling mechanism had good horizontal coupling with the ground, 
but poor vertical coupling. 



NPS Softball Field: 




[ZZ]/ 




Figure 4.4. Field Source Testing Experimental Layout 

At a distance of 1 m, the surface waves generated by the source motion were 
clearly visible, but at 5 m could not be separated from the background noise. This was 
clearly unacceptable, so to increase the chance of detection, improvement in the signal to 
noise ratio was required. First attempts were made by increasing the duty cycle, from 5 
cycles to 20 cycles in increments of five. The increase in signal duration was clearly 
visible at 1 m, but at 5 m the signal was still buried within the background noise. The 
next effort was to increase the drive voltages of the two Bass Shakers. Starting at 5 
Vrms, the drive voltage was increased in increments of five, until the current limit of the 
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power amplifiers was reached at approximately 39 Vrms. Even at 39 Vrms, the surface 
wave motion generated by the source was barely discemable at 5 m. 

As a result of the first field experiment the following conclusions were drawn. 
First, to generate Rayleigh waves of sufficient amplitude to be detected at ranges of at 
least 1 0 m would require that the Bass Shakers be driven at the upper end of their linear 
response voltage limit. Second, that the ground coupling mechanism needed to be 
redesigned so as to enhance the vertical coupling. Furthermore, since there were so many 
channels of data that needed to be observed and ultimately processed, a method of 
recording all the field data was also required. 

D. FIELD TEST SOURCE - MOD II 

Using the lessons learned from the first field experiment, the next major 
milestones were to: 

• Find the operating limits of the Bass Shaker. 

• Modify the ground coupling mechanism to increase the vertical coupling. 

• Incorporate into the field test equipment a method of obtaining multiple data 
channels. 

1. Vibration Source 

The manufacturer specifications that came with the Bass Shakers did not indicate 
any electrical operating limits. The installation guide (Appendix B3) did provide an 
advised maximum continuous power rating of 25 W per shaker. With an input 
impedance of 4 Q and a 25 W power rating, the maximum continuous drive voltage is 
only 10 Vrms. These ratings apply only to its intended use as a car stereo component and 
not for use as a vibration source for a seismic transducer. Consultation with the 
manufacturer provided no further information. 

Without any formal documentation available, the only way to adequately assess 
the operating limits of the Bass Shakers was to see under what operating conditions 
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would failure occur. Sets of failure tests were conducted using the test equipment setup 
as described in Appendix D4. The goal of the failure testing was to estimate the amount 
of signal distortion as a function of frequency and drive voltage. 

The HP35665A Signal Analyzer was used to simultaneously monitor the time 
series and total harmonic distortion of the Bass Shaker using the output signal of a uni- 
axial accelerometer mounted on the shaker case (See Appendix B4 for accelerometer 
specifications). The time series of the accelerometer output was used to visually detect 
any signs of failure or non-linearity as evidenced by changes in signal shape and/or 
duration. An estimate of the signal distortion over the first ten harmonics of the 
fundamental frequency, i.e. the drive frequency, was obtained using a built-in machine 
function. 

The first set of failure tests were conducted by driving the shaker with 20 cycles 
of a sine wave, with a 2.00 sec repetition period, for fundamental frequencies of 100, 
150, 200, 500, and 1000 Hz. While holding the frequency fixed, drive voltages of 10.1, 
20.1, 30.1, and 38.6 Vrms were used, with 38.6 Vrms being the maximum voltage 
attainable without exceeding the current limit of the power amplifier. Figure 4.5 is a 
graphical representation of the signal distortion as a function of frequency and drive 
voltage for a 2.00 sec repetition period. Even though the overall signal distortion was 
low, the accelerometer time series for the 1000 Hz drive frequency showed evidence of 
non-linear behavior. The leading and trailing edges of the accelerometer signal were not 
as detailed as compared to the output signal from the function generator. There was also 
evidence of pulse elongation in the signal. The non-linear behavior was attributed to 
internal heating effects, as evidenced by noticeable case heating observed, but the shaker 
did not fail. 

To stress the Bass Shaker further toward failure required that a more powerful 
amplifier be utilized. To hopefully alleviate the distortion caused by internal heating the 
repetition period was increased to 5.00 sec. A second set of failure tests were conducted 
with the new repetition period using a higher rated amplifier. During this second round 
of testing, the test shaker failed. The shaker was disassembled to locate the cause of the 
failure, but the internal design of the Bass Shaker only made the explanation speculative. 
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Figure 4.5. Signal Distortion for 2.00 sec Repetition Period 



Since the cause of the failure of the first Bass Shaker was unknown, a second 
shaker was setup to complete the evaluation with the new amplifier. Figure 4.6 displays 
the results of the second set of tests. 
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Figure 4.6. Signal Distortion for 5.00 sec Repetition Period 
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Using a repetition period of 5.00 sec seemed to effectively reduce the amount of 
internal heating and resulting distortion. As can be seen from Figure 4.6, the optimal 
drive frequency with respect to minimal signal distortion was 250 Hz. Again since the 
current limit of the power amplifier was exceeded before the shaker failed, it is safe to 
assume that with the current test equipment setup, shaker failure due to input drive 
voltage is not a major concern. 

2. Ground Coupling Mechanism 

To enhance the vertical coupling of the source to ground coupling mechanism, 
two end attachments were mated to the original base as shown in Figure 4.7. The 
attachments were constructed to hold three machine screws each, making the same 45- 




Figure 4.7. Picture of Mod II Base Modifications 
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degree angle with the horizontal surface plane as the Bass Shakers. It was hoped that 
these modifications to the base of the ground coupling mechanism would increase both 
the horizontal and vertical source coupling. 

The upper support brace for the Bass Shakers was also modified. First, it was 
reduced in size to position the shakers closer to the desired 45-degree angle. Second, 
new BNC mounts were incorporated, making the Mod II source more rugged. Figure 4.8 
is a picture of the redesigned Rayleigh wave source. Appendix C3 contains the technical 
drawings for the modifications. 




Figure 4.8. Picture of Field Test Source - Mod II 

Addition of the two end attachments meant that emplacement of the source would 
first require placing the base into the ground in a similar fashion as the Mod I source, 
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then by using a cordless drill boring the six remaining angled machine screws into the 
ground. 



3. Experimental Results 

To resolve the data gathering issue, a SD390 Dynamic Signal Analyzer, made by 
Signal Processing Systems, was added to the test equipment setup described in Appendix 
D3. The SD390 is capable of receiving, displaying, and storing eight channels of 
information simultaneously, in real time. The eight channels provided data gathering 
capabilities for the: 

• Radial and vertical components of the base-mounted accelerometer. 

• Longitudinal, transverse, and vertical components of the two tri-axial 
seismometers. 

Field testing of the Mod II source was also conducted at the softball field, using 
the equipment setup as shown in Appendix D5, but the source experimentation did not 
progress much farther than equipment setup. It was clear from the output signals of the 
accelerometer that the source was not exciting horizontal motion equally as well as it was 
exciting vertical motion. Thus generating the required elliptical motion to selectively 
stimulate Rayleigh waves would be nearly impossible. Attempts to replace the source 
were made, but did not significantly improve the operating conditions. 

The results of the second field experiment prompted the need for further 
refinements to the design of the ground coupling mechanism and to rethink the utility of 
the softball field as the test location. 

E. FIELD TEST SOURCE - MOD III 

It is, of course, important to have good horizontal and vertical coupling to the 
ground for Rayleigh wave stimulation. Generating elliptical particle motion at the source 
requires moving the medium through a combination of both radial and vertical 
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excitations, such that, the semi-major axis of the ellipse is normal to the surface boundary 
and the semi-minor axis parallel to the interface in the direction of propagation. 
Therefore, having the correct design for the ground coupling mechanism is a major 
portion of discrete-mode source development. 



1. Ground Coupling Mechanism 

From the second field experiment it was noted that boring the machine screws 
into the ground with the cordless drill seemed to lock the source in place. The problems 
with the ineffective coupling arose from having to dig holes for the pre-positioned 
vertical machine screws. Hence, it only seemed logical that the source base be modified, 
so that all the screws could be positioned by screwing them into the ground. The Mod III 
base as shown in Figure 4.9 is the result. 




Figure 4.9. Picture of Mod III Base Modifications 
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The Mod III base modifications consisted of the addition of the two larger end 
attachments to the original base and the switch from 3 Vi inch machine screws to 6 inch 
machine screws. Reuse of the original base only served to reduce the turn-around time 
between field tests, since a whole new base did not have to be machined. 

2. Experimental Results 

Due to the equipment difficulties that arose from using the softball field as the test 
location, it was decided that any following ground phase experimentation be conducted 
from the lab in Spanagel Hall. Appendix D6 contains the new test location layout. The 
same equipment setup as described in Appendix D5 was used. 

With the necessary modifications to the ground coupling mechanism completed, 
the next field experiment was used to ensure that the Mod III base provided the required 
horizontal and vertical coupling for Rayleigh wave stimulation. This phase of the source 
testing consisted of verifying, that with the correct phase settings for each of the Bass 
Shakers, both radial and vertical source motion could be individually stimulated. The 
stiffness of the ground to the source motion was quite noticeable, it was only through the 
manual manipulation of both phase and amplitude, that individual excitation of radial and 
vertical motion was possible. 

The fourth, and final, ground phase field experiment was used to investigate the 
axial surface wave field effects for various drive voltages and relative phases of the Bass 
Shakers. Refer to Appendices D5 and D6 for the corresponding equipment and 
experimental setups. The effects to the surface wave field were characterized by 
measuring the coherence between the radial and vertical signal components at the furthest 
seismometer. Coherence on the order of unity would indicate a high degree of reliability 
in the magnitude and phase information of the received signals. The coherence was 
determined using a built-in processing option of the SD390 Signal Analyzer, using the 
average of ten received signals from a 250 Hz, 20 cycle sine wave burst, with a 5.00 sec 
repetition period. Driving voltages of 20. 40, 60, and 86 Vrms and relative phases of 0. 
30, 60, 90, 120, 150, 180 degrees were chosen as the test values. 
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The quantitative results of the experiment are shown in Figure 4.10. As can be 
seen from the figure, the coherence between the radial and vertical components of the 
received signal increases as drive voltage increases, being nearly unity over most phase 
angles for a drive voltage of 86 Vrms. This final ground phase experiment illustrates 
that, the greater the input energy into the medium, the better the stimulation of reliable, 
repeatable, surface wave fields. The reason for the low coherence at 120 degrees is 
unknown. 
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Figure 4.10. Coherence Measurements as a Function of Voltage and Phase 



The time and effort spent in conducting the Ground Phase research provided a 
great amount of detail about the behavior of the field test sources and the associated 
surface wave fields. The choice of the Bass Shakers as the vibration source proved to be 
quite effective in generating the kinds of motion desired for discrete stimulation of 
Rayleigh waves. 
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Since the development of the seismo-acoustic sonar concept, the ultimate goal of 
the combined research efforts was to develop a mine detection system that was capable of 
locating buried ordnance, both in the surf zone and on land. Clearly, the surf zone is the 
more complicated environment to operate in, with regards to both the physics and the 
equipment required. The next chapter describes the development and experimental 
research conducted in this new environment. This portion of the research was identified 
as the Beach Phase. 
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V. DISCRETE-MODE SOURCE DEVELOPMENT AND TESTING 

(BEACH PHASE) 



This chapter covers the beach and surf zone development and experimental 
testing of the discrete-mode, seismic wave source. Included within the chapter are 
individual sections documenting the different phases of the experimentation. The first 
section briefly addresses the final modifications made to the source-to-ground coupling 
mechanism. The next section covers the ambient background noise analysis and ground 
truthing experiments conducted at the test site location. The remaining sections describe 
the research efforts made to assess the ability of the test source to selectively excite 
Rayleigh waves. 

During the ground phase experimentation, research efforts were mainly directed 
at investigating the operating characteristics of the field test sources, in order to find the 
most effective source-to-ground coupling design. Secondly, it was also concerned with 
qualitatively assessing the effects of various source parameters on the axial surface wave 
field. The ground phase research provided the necessary background information that 
was needed to successfully interpret the beach phase experimental results. Progressing 
onto beach and surf zone field experimentation was a major milestone in furthering the 
development of a discrete mode seismo-acoustic source and towards a greater 
understanding of the complex physics involved. 

Using the lessons learned from the ground phase experimentation with respect to 
equipment transportation and logistics, it was evident that working out the bed of a 
pickup truck would not suffice for the beach phase experimentation. Thus, it was 
necessary to find a research platform that was large enough to transport all the equipment 
to the test site and that could provide the necessary electrical power. The mobile home 
owned by the Oceanography Department was the logical and most readily available 
solution. 
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A 



FIELD TEST SOURCE - MOD IV 



From the results obtained during the last set of ground phase experiments, the 
Mod III source-to-ground coupling mechanism, as shown in Figure 4.9, seemed to 
provide the necessary coupling in both the horizontal and vertical directions. As a result 
of its successful testing it was decided that this design be used as a model for the 
fabrication of the final ground coupling mechanism to be used for the beach phase 
experimentation. Figure 5.1 is a picture of the Mod IV ground coupling mechanism. 
Appendix C4 contains the technical drawings for the Mod IV base. 




Figure 5.1. Picture of Mod IV Ground Coupling Mechanism 



46 



To help protect the sources from the harsh beach environment, all joints and 
access holes were filled with silicone. Furthermore, a waterproof accelerometer was also 
added (See Appendix B5 for manufacture specifications). 

B. BACKGROUND NOISE AND GROUND TRUTHING EXPERIMENTS 

Because of the highly complicated nature of the surface wave field, it was 
necessary to obtain a sufficient knowledge base of the ambient acoustic properties of the 
test medium. To obtain the required information meant that ambient background noise 
analyses and ground truthing experiments be conducted at the proposed source testing 
site. The test site chosen was a stretch of beach on the Monterey Bay coastline located 
behind the NPS Beach Research Facilities. The section of beach front area used, as 
shown in Figure 5.2, was approximately 80 feet wide from the surf zone to the sand 
dunes at low tide with a longitudinal extent far exceeding that required for the source 
testing. 
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1 . 



Background Noise Analysis 



In order to make source testing and data interpretation as simple as possible it was 
desired that the source be operated in a frequency range for which competition with the 
background noise would be minimal. Background noise samples were taken at the test 
site using the seismometers as the receivers. Spectral analysis of the noise samples taken 
showed that the strongest noise components resided between 5-20 Hz as shown in Figure 
5.3. Since these frequencies were well below the envisioned operating frequencies for 
the Bass Shakers filtering out the undesired background noise from the input signals 
would be relatively simple. Appendix F contains the MATLAB programs used to 
conduct the spectral analysis. 




2. Ground Truthing 

To obtain an accurate picture of the ambient acoustic properties of the test site a 
series of ground truthing experiments were conducted to estimate the group velocity of 
the Rayleigh wave. As was shown in section II C, Figure 2.6, a vertical point impulse 
generates a wave system comprised of P-, S- and R-waves, as well as other wave types 
not shown. 

To simulate a vertical impulsive source a bowling ball with a self-contained uni- 
axial accelerometer was dropped from an average height of 45 inches, as shown in Figure 
5.4. A bowling ball was used because it is spherically symmetric and impacts the ground 
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the same way on each drop. When the bowling ball impacts the ground several 
radially outward propagating wave systems are created. The impact of the bowling ball 
is sensed by the accelerometer, which triggers a multi-channel data acquisition system to 
record the response of a pre-positioned seismometer. Appendix B6 contains the 
manufacturer specifications for the uni-axial accelerometer and Appendix D7 the wiring 
diagram for the equipment setup. 




Figure 5.4. Picture of Initial Ground Truthing Experiment 

Since a seismometer is a velocity-sensing instrument, the particle motion induced 
by the propagating wave system is converted to an output voltage, which is proportional 
to the particle velocity, which can be used to obtain the wave velocities. The particle 
velocity data obtained from the seismometers were processed using a time versus 
frequency decomposition method. This method involves dividing the seismometer data 
records into equal length blocks with a pre-defined overlap. The midpoint of each block 
can be associated with an arrival time based on the separation between the bowling ball 
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and the seismometer, which in turn can be associated with a wave velocity. The power 
spectral density of each block is calculated and stored in a matrix, which associates 
frequency with velocity. The results of the data processing are then displayed in the form 
of a Gabor plot, named after its inventor Professor Gabor of the University of Delft, the 
Netherlands. The Gabor plot depicts the decomposition of the wave field as a contour 
plot of slowness (inverse speed) as a function of frequency. From the Gabor plots, 
estimates of the propagating velocity of various wave modes can be made. 

Furthermore, since the particle displacement of the beach sediments is 
proportional to the particle velocity, the received seismometer signals can also be used to 
graphically display the average particle motion observed at the seismometer. The 
resulting particle motion can be displayed in the form of a Hankel plot or hodogram. By 
plotting the longitudinal component of motion versus the vertical component as a 
function of time, which is known as a Hankel plot, the motion can be used to determine 
the wave type, i.e. P, S, or R. A hodogram is the two-dimensional form of a Hankel plot. 
Appendix G contains the MATLAB programs used to analyze the ground truthing data. 

The first set of ground truthing measurements were conducted as shown in Layout 
I of Appendix D7. Two series of bowling ball drops were made. The first series was 
conducted approximately 45 feet from the surf zone, in dry sand. The bowling ball was 
dropped at distances of 6, 12, 18, and 24 feet from the seismometer. The second series 
of drops were conducted at the edge of the surf zone using distances of 6, 12, 18, 24, 36, 
and 48 feet. 

During the analysis of the seismometer data it became evident that even a source- 
receiver separation of 48 feet did not allow for the desired time separation of the 
propagating modes in the wavetrain. Thus, an accurate decomposition of the different 
wave types could not be accomplished. Figures 5.5 and 5.6 show the time series of the 
longitudinal and vertical velocity signals received at the seismometer for distances of 24 
feet in dry sand and 48 feet in wet sand, respectively. The time series depicted have been 
adjusted so that 0 sec corresponds with the impact of the bowling ball. Figures 5.7 and 
5.8 are Hankel plots of the particle motion at the seismometer for the associated time 
series. 
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Recalling that Rayleigh waves induce elliptical particle motion, its clear from the 
Hankel plots that the motion shown at these distances is due to the presence of more than 
just the Rayleigh wave. Thus to adequately estimate the various wave speeds would 
require allowing greater time for wave separation to occur. This could only be 
accomplished by using larger separation distances between the seismic source and the 
seismometer. 

To promote wave separation, the same experiment was conducted a second time 
using much larger distances. Figure 5.9 and Layout II in Appendix D7 show the 
experimental layout for the second phase of ground truthing experiments. Separation 
distances of 50, 70, 90, 1 10, and 130 feet were selected. Also, a second seismometer was 
placed at a distance of 10 feet behind the first to allow for the collection of five additional 
data sets at intermediate distances. Furthermore, it was decided that the ground truthing 
efforts be concentrated in the drier areas of the beach since the Bass Shakers were not 
waterproof. 



Longitudinal Time Series of Geophone Signal for File: dry24ft1 





Figure 5.5. Time Series of Seismometer Velocity Signals for a 24 ft. 
Source-Receiver Separation Distance 
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x 10 ~ 3 Longitudinal Time Series of Geophone Signal for File: wet48ft1 





Figure 5.6. Time Series of Seismometer Velocity Signals for a 48 ft. 
Source-Receiver Separation Distance 




Figure 5.7. Hankel Plot of Seismometer Velocity Signals for a 24 ft. 
Source-Receiver Separation Distance 
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Figure 5.8. Hankel Plot of Seismometer Velocity Signals for a 48 ft. 
Source-Receiver Separation Distance 




Figure 5.9. Picture of Second Ground Truthing Experiment 
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Increasing the separation distance between the drop location and the first 
seismometer greatly increased the ability to decompose the wave field. It was not until 
the data files for distances of 70 feet and greater were analyzed that modal separation in 
the Gabor plots became clear. Figure 5.10 is the Gabor plot for a separation distance of 
70 ft (21.3 m). Figure 5.10 shows the existence of at least two wave modes, the slower 
mode being the Rayleigh wave as determined by analysis of the hodogram of the time 
series velocity signals. The group velocity of the Rayleigh wave at 40 Hz is 
approximately 88 m/sec (0.0114 sec/m) and the second mode is another wave with a 
group velocity of approximately 1 14 m/sec (0.0088 sec/m) at 60 Hz.. The identity of the 
second wave is unknown since its speed is not consistent with that expected for a shear 
wave in this sediment. 
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Figure 5.10. Gabor Plot for a 70 ft. Source-Receiver Separation Distance 
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The forces generated by the bowling ball impact with the beach sediments seems 
to generate larger amplitude body waves, as compared to the point source discussed in 
section II C. This is clearly evident in the time series for the seismometer longitudinal 
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section II C. This is clearly evident in the time series for the seismometer longitudinal 
velocity signal displayed in Figure 5.1 1. Even at a distance of 130 feet, the longitudinal 
particle motion caused by the body waves is comparable to that of the Rayleigh wave. 
The displacements in the vertical direction are almost entirely dominated by the Rayleigh 
wave. 

Figures 5.12 and 5.13 are the Hankel plot and hodogram, respectively, for the 
time series of Figure 5.1 1. These figures serve to illustrate how much more prevalent is 
the retrograde elliptical motion of the Rayleigh wave at these larger distances. Figure 
5.14 is the Gabor plot of the wave decomposition at a range of 130 ft (39.6 m). As can be 
seen from the figure, two propagating wave modes are present in the velocity range 60- 
250 m/sec. From the Gabor plot, the Rayleigh wave speed is estimated at approximately 
89 m/sec (0.01 12 sec/m) and the second wave has a speed of approximately 109 m/sec 
(0.0092 sec/m). The ratio of the Rayleigh and shear speeds obtained from this plot are 
within 10% of the theoretical ratio predicted by the theory of elasticity. 



Longitudinal Time Series of Geophone Signal for File: bl 301402 





Figure 5.1 1. Time Series of Seismometer Velocity Signals for a 130 ft. 
Source-Receiver Separation Distance 
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Figure 5.12. Hankel Plot of Seismometer Velocity Signals for a 130 ft. 
Source-Receiver Separation Distance 




Figure 5.13. Hodogram of Seismometer Velocity Signals for a 130 ft. 
Source-Receiver Separation Distance 
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Figure 5.15 shows an aggregate display of the accelerometer trigger signal (lowest 
trace) and the seismometer vertical velocity versus time recorded for the bowling ball 




Figure 5.15. Rayleigh Wave Progression as a Function of Space and Time 
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drops between 50 ft (15.2 m) and 130 ft (39.6 m) range. The behavior of the velocity 
amplitude with range should not be considered indicative of what would be observed for 
an end-on array of five seismometers receiving a single signal, as site conditions probably 
changed over the duration of the data collection. The large vertical amplitude of the 
Rayleigh wave is easily distinguishable in the more distant traces. The extracted group 
velocity is approximately 90 m/sec (see straight line in Figure 5.15), which agrees well 
with speeds estimated from the Gabor plots. 



C. DISCRETE-MODE EXCITATION 

With the necessary ground truthing experiments completed, the next phase of the 
research was to begin testing the Mod IV field test source (See Figure 5.16). The goal of 



Figure 5.16. Picture of Mod IV Field Test Source 
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the research conducted from this point on was to prove that by inducing the appropriate 
elliptical particle motion into the beach sediments at the source, Rayleigh waves might be 
preferentially excited. Based on the theory for surface waves in a homogeneous, 
isotropic, elastic half-space, the particle motion induced by Rayleigh waves at the 
boundary surface is retrograde elliptical, with the vertical motion approximately l -'A 
times the longitudinal motion. 

1. Propagation Characteristics 

Before moving onto accessing the Mod IV source’s ability to selectively excite 
Rayleigh waves, some brief experimentation was conducted to see if the beach site itself 
had a preferential direction for energy propagation. To investigate this, the Mod IV field 
test source was positioned with two seismometers at equal distances, in the longitudinal 
direction, on both sides of the source, as shown in Figure 5.17 and layout I of Appendix 




Figure 5.17. Picture of Wave Propagation Experiment 
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D8 (Also see Appendix D8 for associated test equipment setup). 

The source as shown in Figures 5.16 and 5.17 was positioned so that, facing 
seaward, the port Bass Shaker faced in southerly direction along the beach and the 
starboard Bass Shaker to the north (See Figure 5.2 for orientation). Analysis of the 
shaker motion indicated that for an equal amplitude, positive input voltage, a negative 
phase for the port shaker and a positive phase for the starboard would generate retrograde 
elliptical source motion with respect to the southern positioned geophone. Reversing the 
phase would then induce a prograde elliptical motion. 

For a homogeneous, isotropic, elastic medium the ratio of the longitudinal to 
vertical particle displacements for the Rayleigh wave at the boundary surface is a well- 
defined function of the Poisson's ratio. But since the Poisson’s ratio for the beach 
sediments was not known, the specific ratio to be used in driving the test source was also 
unknown. As a result, it was decided that initial source testing be done using a 90 degree 
phase difference between the two shakers. This was accomplished by setting the port 
shaker to ^45 degrees and the starboard shaker to +45 degrees, or vice versa. To 
distinguish between the two cases the following convention will be used: a -45 degree 
phase for the port shaker will be denoted by —45P and a +45 degree phase by +45P. 

Since the Bass Shakers suffered the least distortion at 250 Hz, (See section 
IV.D.l, Figure 4.6), this frequency was chosen as the initial drive frequency. A 50 V 
peak, 20 cycle, 250 Hz, narrow-band filtered tone burst was used to drive the shakers. 
Even though this frequency was optimal for the shakers, it became immediately evident 
during data collection that this drive frequency was not optimal for energy propagation 
along the beach. As a result, it was decided to reduce the drive frequency to 100 Hz. It 
has been measures that, for this frequency, the harmonic distortion is less than 1% for a 
50 Vrms drive voltage. Figures 5.18 and 5.19 are the Hankel plots for the seismometer 
signals at a range of 30 ft using a 50 V peak, 20 cycle, 100 Hz, narrow-band filtered tone 
burst and a -45P phase configuration. Figure 5.18 depicts the particle motion as 
observed at the southern seismometer and Figure 5.19 for the northern seismometer. In 
both figures the elliptical particle motion is prograde with respect to the seismic source 
and clearly shows the presence of the Rayleigh wave, but the northerly propagating wave 
is much more uniform. It can also see from the figures that the amplitude of the northerly 
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propagating wave is larger by a factor of two. Appendix H contains the MATLAB 
programs used to analyze the data gathered during the beach phase experimentation with 
the Mod IV source and to generate the accompanying figures. 




time fsecl 

Figure 5.18. Hankel Plot of Southern Seismometer Velocity Signals for a 30 ft. 
Source-Receiver Separation Distance (-45 P) 
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Figure 5.19. Hankel Plot of Northern Seismometer Velocity Signals for a 30 ft. 
Source-Receiver Separation Distance (-45P) 
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Figures 5.20 and 5.21 are the Hankel plots for the same seismometer orientation 
and range, but for a +45P phase setting, which also produced an elliptical particle motion 
that was prograde with respect to the field test source. Again we see that the Hankel plot 
for the southern seismometer (Figure 5.20) is non-uniform and has a smaller amplitude as 
compared to the plot for the northern seismometer (Figure 5.21). Based on these figures, 
it would seem that the beach sediments either have a preferential direction of energy 
propagation or that the force generated by the motion of test source is different in the two 
directions. 




Figure 5.20. Hankel Plot of Southern Seismometer Velocity Signals for a 30 ft. 
Source-Receiver Separation Distance (+45P) 
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Figure 5.21 . Hankel Plot of Northern Seismometer Velocity Signals for a 30 ft. 

Source-Receiver Separation Distance (+45 P) 

The next set of the beach experiments was used to investigate the typical 
propagation distances that could be achieved with the discrete mode source. Since the 
results of the previous source tests were not available at the time of conducting this phase 
of the beach experimentation, the choice for follow-on seismometer position was 
completely arbitrary. 

The two seismometers were initially positioned at distances of 30 and 45 ft and 
leapfrogged out to final distances of 90 and 105 ft. as shown in Layout II of Appendix 
D8. Using the same 50 V peak, 20 cycle, 100 Hz tone burst, sufficient Rayleigh wave 
signal strength was observed out to a distance of 90 ft. which far exceeded any of our 
initial expectations. Propagation distances of this order are just one of the features that 
make seismic sonar research so extremely important to mine warfare applications, since 
the maximum propagating distances for the best ground penetrating radar and 
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magnetometers is only 25 ft. Figure 5.22 shows an aggregate display of the source 
vertical acceleration (lowest trace) and the seismometer vertical velocity versus time at 
ranges from 30 ft (9.14 m) to 90 ft (27.43 m). Each seismometer signal trace has been 
normalized by the source vertical acceleration for that trace, as measured by the base 
mounted accelerometer. The behavior of the velocity amplitude with range should not be 
considered indicative of what would be observed for an end-on array of five 
seismometers receiving a single signal, as site conditions (e.g. tide, surface moisture, etc.) 
were evolving during the course of the data collection. The Rayleigh wave is clearly 
visible as the large amplitude wave train in each trace. The extracted group velocity is 
approximately 86 m/sec (see straight line in Figure 5.22), which agrees very well with 
that observed during ground truthing measurements. 
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Figure 5.22. Rayleigh Wave Progression as a Function of Space and Time 

Another interesting observation from this phase of the beach experimentation was 
that analysis of the particle motion observed at the seismometers showed prograde 
elliptical motion with respect to the source for all ranges, seismometer orientations, and 
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phase configurations. This result was unexpected, but can be explained. Assuming the 
beachfront to be a nearly homogeneous, isotropic, elastic medium, then the Rayleigh 
wave speed is independent of frequency (Refer to section II B). For a Rayleigh wave 
speed of 86 m/sec, a 100 Hz drive frequency corresponds to a Rayleigh wavelength of 
0.86 m. Experimentation shows that the shift from retrograde to prograde elliptical 
motion for sediments occurs at approximately 0.1 X r or 8.6 cm (3.4 in) [Ref. 9]. Since the 
seismometers nominally sense the average particle motion and are nearly 25.4 cm (10 in) 
in length, most of their bulk resides below this depth. Thus, the observed motion is 
consistent with theory. 

In order to verify this observation, another set of experiments were conducted 
using a 50 Hz drive frequency. At 50 Hz, assuming the same Rayleigh wave speed of 86 
m/sec, the Rayleigh wavelength is 0.34 m. For this wavelength the corresponding turning 
point occurs at a depth of 17.2 cm (6.8 inches). To help eliminate any chance of 
misinterpreting the data the seismometers were only buried to a depth of 12.7 cm (5 in). 
Figures 5.23 and 5.24 are the Hankel plots of the seismometer velocity signals at 60 feet 
for drive frequencies of 1 00 Hz and 50 Hz, respectively. 

As can be seen from the figures the sense of the particle motion is clearly 
different. Analysis of the Hankel plots showed that the particle motion as shown in 
Figure 5.23 was prograde elliptical and that of Figure 5.24 was retrograde elliptical. This 
illustration serves to show that accounting for the vector nature of the Rayleigh wave is 
extremely important, not only for discrete-mode propagation, but also for selective 
reception. 
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Figure 5.23. Hankel Plot of Seismometer Velocity Signals for a 60 ft. 
Source-Receiver Separation Distance (100 Hz) 
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Figure 5.24. Hankel Plot of Seismometer Velocity Signals for a 60 ft. 
Source-Receiver Separation Distance (50 Hz) 
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2 . 



Selective Excitation 



One of the goals of this thesis was to attempt to estimate how much of the input 
mechanical energy was being converted into Rayleigh waves. But because of the 
unpredictable nature of the source motion, this was not possible. It was assumed with the 
source design chosen and the information obtained on the operating characteristics of the 
Bass Shakers that it would be relatively easy to generate elliptical source motion. During 
the beach phase experimentation it was discovered that it was nearly impossible just to 
create circular motion at the source. Using equal amplitude drive voltages for the shakers 
and a — 45P phase configuration should have produced retrograde circular motion and a 
+45P, prograde circular motion. But, as can be seen from the hodograms of the base 
mounted accelerometer for the two cases. Figures 5.25 and 5.26, this was clearly not the 
result. 

0.1 F 1 I T 1 I =1 




-o.i - 

i i i i i 

- 0.1 - 0.05 0 0.05 0.1 

longitudinal amplitude [V] 



Figure 5.25. Hodogram of Accelerometer Signal (100 Hz, -45P) 
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Figure 5.26. Hodogram of Accelerometer Signal (100 Hz, +45P) 



Analysis of the hodograms from all the data files showed that the motion 
produced at the source depended explicitly on voltage, frequency, and phase. 
Furthermore, there is also a strong dependence on the elastic properties of the 
propagating medium. Examination of the accelerometer hodograms showed that 
longitudinal stimulation of the beach sediments was much more effective than vertical 
stimulation. 

Follow-on experiments are planned in which it will be attempted to produce more 
optimum source motion. Even though the source motion was clearly not optimum, this 
did not seem to play a major part in stimulating the desired Rayleigh wave, at least at 100 
Hz. Figure 5.27 is the Hankel plot of the seismometer velocity signals at 60 ft for a 20 V 
peak, 20 cycle, 100 Hz, tone burst with no phase difference between the shakers. Figure 
5.28 displays the portion of the time record corresponding to the Rayleigh wave and the 
ratio of the longitudinal to vertical amplitude over the same time period. 
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Figure 5.27. Hankel Plot of Seismometer Velocity Signals for a 60 ft. 
Source-Receiver Separation Distance 
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Figure 5.28. Ratio of Longitudinal and Vertical Rayleigh Wave Components 
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The uniformity of the wave is remarkable. Similar calculations were done for the 
seismometer velocity signals at 60 feet when phase differences of — 45P, +45P, -32P, and 
+32P were used. In all four cases, similar looking Rayleigh waves were excited and the 
ratio of the wave components was nearly identical. 

From analysis of the data files gathered during the beach phase experimentation it 
was obvious that the Rayleigh wave was the most significant disturbance along the 
interface boundary. At distances of 45 ft and smaller it appeared that no matter what 
phase configuration was used, all wave types were present, with the Rayleigh wave and 
vertically polarized shear wave the most dominant. But since body waves along the 

surface decay like y and Rayleigh waves decay like , the influence of the body waves 
was negligible at distances over 60 ft. 
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VI. 



CONCLUSIONS AND RECOMMENDATIONS 



This final chapter briefly synopsizes the noteworthy observations made during the 
conduct of the research described in this thesis and the conclusions that can be drawn 
from these observations. It also includes recommendations for follow-on research. 

The major goal of this research project was to evaluate the concept of a discrete- 
mode interface (Rayleigh or Scholte) wave source. By selectively exciting the desired 
interface waves, it was theorized that significant reverberation suppression could be 
achieved, thus increasing the signal-to-noise ratio of a target echo in a seismic sonar 
detection system. 

Based on the precepts from the theory of elasticity and the known characteristics 
of the desired interface waves, two different source concepts were designed and tested. 
The first discrete-mode source, i.e., the concept source as shown in Figure 4.1, was used 
to qualitatively probe in a laboratory environment the surface wave field created under 
various source motions. Experimentation with the concept source showed that a complex 
relationship existed between input source parameters, such as, frequency, phase, and 
driving voltage, and the propagating characteristics of the medium, and the resulting 
surface wave field. The limitations of the concept source and those imposed by the test 
tank made development of a field test source a necessity if any real progress was to be 
made. 

To design a source that could be used under varying field conditions first required 
identifying a suitable vibration source and then constructing a matching source-to-ground 
coupling mechanism. Figure 4.3 shows the first version of the field test source. Again 
the beginning phases of the field experimentation were directed at evaluating the effects 
of varying source motions on the surface wave field. Results obtained during the various 
field tests lead to several modifications of the source-to-ground coupling mechanism until 
arriving at what was thought to be the best design (See Figure 5.16). It was with this 
final source design that all the beach phase testing was conducted. 

Based on the known features of Rayleigh waves, it was believed that by 
stimulating the sediments at the source in an elliptical motion, Rayleigh waves could be 
selectively excited. To test this theory a series of experiments were conducted in which 



71 



the field test source was driven under several different modes of excitation. During these 
experiments it was observed that multiple types of interface waves and body waves were 
generated no matter what input source parameters were used. At source-to-receiver 
ranges greater than approximately 1 5 m (45 ft), the Rayleigh wave was clearly the largest 
interface disturbance, while the effects of the undesired waves quickly became negligible 
and were virtually undetectable at ranges greater than approximately 25 m (75 ft). 
Evaluation of the base-mounted accelerometer on the source showed that the field test 
source did not produce the elliptical motion expected for the applied input parameters. 
Because of this, it is difficult to determine whether the excitation of the undesired waves 
was due to the non-optimum nature of the source or that these undesired waves will 
always be present in practice. However, it was found that the medium itself acted as a 
selective filter for the interface waves after a few tens of meters of propagation. 

Based on the experimental results obtained there are two major areas in which 
future discrete-mode excitation research efforts should be concentrated. The first area 
deals with correcting the deficiencies noted with the Bass Shakers and the source design. 
The Bass Shakers were not designed to be operated in the fashion utilized during the 
source testing. Thus, higher quality, greater power-handling exciters (e.g. shakers, 
actuators) must be incorporated into a new source. New sources must be identified that 
can be driven harder, but experience less distortion. Furthermore, complex source 
geometries such as the 45-degree orientation of the shakers used in the field test sources 
should be avoided. Such configurations simultaneously generate both longitudinal and 
vertical seismic motion. Future source designs should utilize simple vertical and 
horizontal excitation. This would serve two purposes, the first being that such an 
orientation still allows for the generation of elliptical motion at the source, and second 
that each component of the motion can be excited individually to study the resulting 
effects. Work has already begun in designing and constructing a new source made of 
orthogonal, servo-controlled, 20 lbf, linear actuators. 

The second area deals with re-examining the experimental setups in order to 
better sample the various seismic wave types. This would serve to increase the 
researcher’s knowledge of the different excitation potentials of the possible wave types 
and their associated excitation depths. Along these lines more experiments are needed to 



72 



further explore source control for optimal wave generation as a function of frequency and 
excitation depth. To better accomplish this task efforts should be made to computerize 
the source control so that the source motion can be monitored on a real-time basis. 

Even though the results of this thesis were not as informative as initially desired 
the basic concept of discrete-mode excitation is promising. An important 
accomplishment of this thesis was to get the research out of the laboratory and onto the 
beach where the follow-on efforts can begin. 
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APPENDIX A. THEORETICAL PROPERTIES OF RAYLEIGH WAVES IN A 
HOMOGENEOUS, ISOTROPIC, ELASTIC HALF-SPACE 



This appendix contains the various MATLAB programs that were used to 
reproduce the theoretical graphs contained within the reference texts on Rayleigh waves, 
some of which appear as figures within the body of the thesis. It also contains other 
programs written and their output, which provide further visual aid into the behavior of 
Rayleigh waves. 
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APPENDIX Al. VARIATIONS IN WAVE SPEED RATIOS AS A FUNCTION OF 

POISSON’S RATIO 



Filename: wavespeeds.m 
Written by: F.E. Gaghan 
% Date Last Modified: 22 December 1997 

% Purpose: This program numerically solves the Rayleigh wave equation 
% in an isotropic elastic half space. The wave equation 

% solutions are expressed as the ratio of the Rayleigh wave 

% speed to the shear wave speed as a function of the 

% Poisson's ratio. The ratio of the three wave types to the 

% shear wave speed is displayed graphically for Poisson's 

% ratios in the interval [0,1/2]. 

% Initialize required storage vectors 
clear all; 

speed_ratiosqd_root 1= [ ] ; speed_ratiosqd_root2= [ ] ; 
speed_ratiosqd_root3= [ ] ; 

longsqd_ratiol= [ ] ; longsqd_rat io2= [ ] ; longsqd_rat io3= [ ] ; 
vert sqd_ratiol= [ ] ; vertsqd_ratio2= [ ] ; vertsqd_ratio3= [ ] ; 
Rayleighsqd_root= [ ] ; 

% For loop used to run through the desired Poisson's ratios 
j=l; 

for nu= [ 0 : 0 . 01 : 0 . 5 ] ; 

gammasqd (j) = (l-(2*nu) ) / (2- (2*nu) ) ; 
inv_gammasqd ( j ) =1 /gammasqd ( j ) ; 

% Calculate the coefficients of the Rayleigh wave equation which is 
% cubic in the wave speed ratio 

p=-8; 

q=24- ( 16 ^gammasqd ( j ) ) ; 

r= ( 16*gammasqd ( j ) ) -16; 

a= ( ( 3*q) - (p A 2 ) ) /3 ; 

b=( (2* (p A 3) )-(9*p*q)+(27*r) )/27; 

f =sqrt ( ( (b A 2 ) / 4 ) + ( (a A 3) /27) ) ; 

g=~ (b/2 ) +f ; 

h=- (b/2 ) -f ; 

% Perform necessary manipulations of the coefficients to acquire 
% values in the correct quadrants when taking cube roots of complex 
% quantities 

if ( f >0 ) & (g<0) 

A= (abs (g) A (1/3) ) *-l; 
else 

A=g A (1/3) ; 

end;%end-if ( f > 0 ) & (g<0) 

if (f>0) & (h<0 ) 

B= ( abs (h) A (1/3) ) *-l; 
else 

B=h A (1/3) ; 

end;%end-if (f>0) & (h<0) 



76 



% Calculate and store the three roots of the square of the wave 
% speed ratio 
xl-A+B; 

x2=(- (A+B) /2) + ( ( (A-B) /2) *sqrt (-3) ) ; 
x3=(- (A+B) / 2 ) — ( ( (A-B) /2) *sqrt (-3) ) ; 
yl=xl- (p/3) ; 
y2=x2- (p/3) ; 
y3=x3- (p/3) ; 

speed_ratiosqd_root 1= [ speed_ratiosqd_root 1 ; yl] ; 
speed_ratiosqd_root2= [ speed_ratiosqd_root2 ; y2] ; 
speed_ratiosqd_root 3= [ speed_ratiosqd_root3 ; y3] ; 

% Calculate and store the ratio of longitudinal attenuation 
% coefficient squared to the Rayleigh wave number squared 
longsqd_rat iol= [ longsqd_ratiol ; 1- ( gammasqd ( j ) *yl ) ] ; 
longsqd_rat io2= [ longsqd_ratio2 ; 1- ( gammasqd ( j ) *y2 ) ] ; 
longsqd_rat io3= [ longsqd_ratio3 ; 1- ( gammasqd ( j ) *y3) ] ; 

% Calculate and store the ratio of vertical attenuation coefficient 
% squared to the Rayleigh wave number squared 
vert sqd_rat iol= [ vertsqd_ratiol ; 1-yl] ; 
vert sqd_rat io2= [ vertsqd_ratio2 ; l-y2] ; 
vertsqd__ratio3= [ vertsqd_ratio3 / l-y3] ; 

j=j+l; 

end;%end-for nu 

% Store the vectors into matrices 

speed_ratiosqd_roots= [speed_ratiosqd_root 1 , speed__rat iosqd_root2 , 
speed_rat iosqd_root 3] ; 

longsqd_rat ios= [ longsqd_ratiol , longsqd_rat io2 , longsqd_rat io3 ] ; 
vert sqd_ratios= [ vertsqd_ratiol , vertsqd_ratio2 , vertsqd_rat io3 ] ; 

% Perform required logical comparisons to keep only the roots that 
% satisfy the Rayleigh wave criteria 

suitable= ( longsqd_ratios>=0 ) & ( vertsqd_rat ios>=0 ) ; 

for m=l : ( j -1) 
for n=l : 3 

if suitable (m, n) 

Rayleighsqd__root= [ Rayleighsqd_root ; speed_ratiosqd_root s (m, n) ] ; 
end;%end-if suitable 
end;%end-for n 
end; %end-for m 

Rayleigh_root=sqrt ( Rayleighsqd_root ) ; 
nu= [0:0.01:0.5] ; 
shear_speed=ones ( 1, length (nu) ) ; 

% Plot the ratio of Rayleigh wave speed to shear wave speed versus 
% Poisson’s ratio 
figure ( 1 ) ; 

plot (nu, shear_speed, nu, Rayleigh_root , nu, inv_gammasqd) ; 
axis ( [ 0 , 0 . 5 , 0 , 7 ] ) ; 
grid on; axis square; 

xlabel (' Poisson' s ratio’); ylabel('wave speed/shear wave speed'); 
text ( . 1 1 , 1 . 2 , ’Shear ') ; text (.11,0.65, ’ Rayleigh ' ) ; 
text (0.305, 3. 2, ' Longitudinal ' ) ; 
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title (’ Variation in Wave Speed Ratios as a Function of Poisson's 
Ratio ’ ) ; 



Variation in Wave Speed Ratios as a Function of Poisson's Ratio 
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APPENDIX A2. VARIATION IN RAYLEIGH WAVE SPEED AS A FUNCTION 

OF POISSON’S RATIO 



% Filename: rayleighspeed . m 
% Written by: F.E. Gaghan 
% Date Last Modified: 02 December 1997 

% Purpose: This program numerically solves the Rayleigh wave equation 
% in an isotropic elastic half space. The wave equation 

% solutions are expressed as the ratio of the Rayleigh wave 

% speed to the shear wave speed as a function of the 

% Poisson's ratio. The wave speed ratio is displayed 

% graphically for Poisson's ratios in the interval [-1,1/2]. 

% Initialize required storage vectors 
clear all; 

speed__ratiosqd__root 1= [ ] ; speed_rat iosqd_root 2= [ ] ; 
speed_ratiosqd__root3= [ ] ; 

longsqd_ratiol= [ ] ; longsqd_rat io2= [ ] ; longsqd__rat io3= [ ] ; 
vertsqd_ratiol= [ ] ; vert sqd_ratio2= [ ] ; vertsqd_rat io3= [ ] ; 
Rayleighsqd_root= [ ] ; 

% For loop used to run through the desired Poisson's ratios 
j=l; 

for nu= [ 0 : 0 . 05 : 0 . 5 ] ; 

gamma sqd= (1- (2*nu) ) / (2- (2*nu) ) ; 

% Calculate the coefficients of the Rayleigh wave equation which is 
% cubic in the wave speed ratio 

p=-8; 

q=24- ( 16*gammasqd) ; 

r=(16*ganunasqd) -16; 

a= ( ( 3*q) - (p~2 ) )/3; 

b=( (2* (p A 3 ) ) - ( 9*p*q) + (27*r) ) /27; 

f =sqrt ( ( (b"2) /4) +( (a^3) /27) ) ; 

g=- (b/2) +f ; 

h=- (b/2) -f ; 

% Perform necessary manipulations of the coefficients to acquire 
% values in the correct quadrants when taking cube roots of complex 
% quantities 

if (f>0) & (g<0) 

A= ( abs (g) A (1/3) ) *-l ; 
else 

A=g A (1/3) ; 

end;%end-if (f>0) & (g<0) 

if (f>0) & ( h< 0 ) 

B= (abs (h) A (1/3) ) * - 1 ; 
else 

B=h A (1/3) ; 

end;%end-if ( f >0 ) & (h<0) 
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% Calculate and store the three roots of the square of the wave 
% speed ratio 
xl=A+B ; 

x2= (- (A+B) /2) + ( ( (A-B) / 2 ) *sqrt (-3) ) ; 
x3= (- (A+B) /2) - ( ( (A-B) / 2 ) *sqrt (-3) ) ; 
yl=xl- (p/3 ) ; 
y2=x2- (p/3) ; 
y3=x3- (p/3 ) ; 

speed_rat iosqd_root 1= [ speed_ratiosqd_root 1 ; yl] ; 
speed_ratiosqd_root2= [speed_ratiosqd_root2 ; y2 ] ; 
speed_ratiosqd_root 3= [speed_ratiosqd_root 3/ y3] / 

% Calculate and store the ratio of longitudinal attenuation 
% coefficient squared to the Rayleigh wave number squared 
longsqd_ratiol= [ longsqd_ratiol ; 1- ( gamma sqd*yl ) ] ; 
longsqd_rat io2= [ longsqd_ratio2 ; 1“ ( gamma s qd * y2 ) ] ; 
longsqd_ratio3= [longsqd_ratio3; 1- (gammasqd*y3) ] ; 

% Calculate and store the ratio of vertical attenuation coefficient 
% squared to the Rayleigh wave number squared 
vertsqd_ratiol= [ vert sqd_ratiol ; 1-yl] ; 
vertsqd_ratio2= [ verts qd_r at io2 ; l-y2] ; 
vertsqd_ratio3= [ vertsqd_ratio3 ; l-y3] ; 

j=j+l; 

end;%end-for nu 

% Store the vectors into matrices 

speed_ratiosqd_roots = [speed_ratiosqd_root 1 , speed_rat iosqd_root2 , 
speed_ratiosqd_root3 ] ; 

longsqd_ratios= [ longs qd_r at iol , longs qd_r at io2 , longsqd_rat io3] ; 
vertsqd_ratios= [ verts qd_r at iol , vert sqd_rat io2 , vertsqd_rat io3] ; 

% Perform required logical comparisons to keep only the roots that 
% satisfy the Rayleigh wave criteria 

suitable= (longsqd_ratios>=0) & ( vert sqd_ratios>=0 ) ; 

for m=l : ( j-1) 
for n=l:3 

if suitable (m, n) 

Rayleighsqd_root= [Rayleighsqd_root ; speed_rat iosqd_roots (m, n) ] ; 
end;%end-if suitable 
end;%end-for n 
end; %end-for m 

Rayleigh_root=sqrt (Rayleighsqd_root ) ; 
nu= [0:0.05:0.5] ; 

% Plot the ratio of Rayleigh wave speed to shear wave speed versus 
% Poisson’s ratio 
figure ( 1 ) ; 

plot (nu, Rayleigh_root ) ; 
grid on; 

xlabel (' Poisson ' s ratio 1 ); ylabel (' Rayleigh wave speed/shear wave 
speed ’ ) ; 

title ( 1 Variation in Rayleigh Wave Speed Ratio as a Function of 
Poisson's Ratio’); 
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Variation in Rayleigh Wave Speed Ratio as a Function of Poisson's Ratio 
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APPENDIX A3. ELLIPTICAL PARTCILE MOTION INDUCED BY RAYLEIGH 

WAVE 



Filename: rayleighellipses . m 

Written by: F.E. Gaghan 

Date Last Modified: 02 December 1997 

Purpose: This program is used to graphically display the varying 

particle motions induced by propagating Rayleigh waves in an 
elastic medium. Specifically it is used to show how the 
particle motion begins as retrograde elliptical and 
progresses to prograde elliptical with depth. 

Initialize required space and wave parameters 
clear all; 
k= (2*pi*100 ) / 3 0 0 ; 
w=k*300; 

t= [0:0. 1:2] /200; 
wt=w*t ; 
wt_t 0=0 ; 

C=0 . 2 ; 

Set special plot commands 

axes ( 'box’ , 'on' ) ; hold on; grid on; axis equal; 

Outer for loop used to step through five various depths along the 

z-axis within the medium 

for x3=[0, -0 . 3, -0 . 6, -0 . 9, -1 . 2] 

% Inner loop used to step through the longitudinal ranges in the 
% direction of propagation along the x-axis 
for xl= [ 0 , 0 . 5 , 1 , 1 . 5, 2, 2 . 5] ; 

% Algebraic manipulations required to calculate the longitudinal 
% particle displacement 

al=exp (0.8475*k*x3)-0. 5773*exp (0.3933*k*x3) ; 
cl=sin (wt- ( k*xl) ) ; 
ul=C* (al*cl) ; 
xlul=ul+xl ; 

% Algebraic manipulations required to calculate the vertical 
% particle displacement 

a3=-0 . 8 4 7 5 *exp (0.8475*k*x3)+1.4 67 9* exp (0.3933*k*x3) ; 
c3=cos ( wt- ( k*xl ) ) ; 
u3=C + (a3 + c3 ) ; 
x3u3=u3+x3 ; 

% Plot elliptical particle motion 
figure ( 1 ) ; 

plot (xlul, x3u3, 1 g . ' ) ; 
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vertical depth [m] 



% Algebraic manipulations required to calculate the wave motion 
% in the longitudinal direction 
xl_tO= [0:0.05:2.5] ; 
last=length (xl_t0) ; 
x3_t0=x3*ones ( 1 , last ) ; 
cl_t0=sin (wt_t0- ( k. *xl_t0) ) ; 
ul_tO=C* (al . *cl_t0) ; 
xlul_t 0=ul_t 0+xl_t0 ; 

% Algebraic manipulations required to calculate the wave motion 
% in the vertical direction 
c3_t0=cos (wt_t0- (k. *xl_t0) ) ; 
u3_tO=C llr (a3 . *c3_t0 ) ; 
x3u3_t0=u3_t 0+x3 ; 
plot (xlul_t0 / x3u3_t0, ' b ' ) ; 

% Plot Rayleigh wave motion 

quiver (xl_t0 (1:10: last) , x3_t0 ( 1 : 10 : last ) , ul_t0 ( 1 : 10 : last ) , ... 
u3_t0 (1:10: last) ,0.125, f k') ; 

end;%end-for xl 
end;%end-for x3 

xlabel (' longitudinal range [m] ’ ) ; ylabel (' vertical depth [m] ' ) ; 
tit le (' Elliptical Particle Motion Induced by Rayleigh Wave’); 



Elliptical Particle Motion Induced by Rayleigh Wave 
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APPENDIX A4. LONGITUDINAL AND VERTICAL PARTICLE 

DISPLACEMENTS 



% Filename: rayleighdisplacement s . m 
% Written by: F.E. Gaghan 
% Date Last Modified: 02 December 1997 

% Purpose: This program is used to graphically display the horizontal 
% and vertical particle displacements induced by a propagating 

% Rayleigh wave for a specified value of Poisson's ratio. The 

% curves displayed are the particle displacement amplitudes 

% normalized by the vertical particle displacement at the 

% surface as a function of depth normalized by the Rayleigh 

% wave number. 

% Initialize required storage vectors 
clear all; 
nu=0 . 25; 

gamma sqd= (1- (2*nu) ) / (2- (2*nu) ) ; 

% Calculate the coefficients required to solve the Rayleigh wave 
% equation which is cubic in the wave speed ratio 

p=-8; 

q=24- ( 16* gamma sqd) ; 

r= ( 16*gammasqd) -16; 

a = ( ( 3 * q ) - (p A 2) ) /3; 

b=( (2* (p A 3) ) - (9*p*q)+ (27 *r) ) / 27 ; 

f =sqrt ( ( (b A 2) /4) + ( (a A 3) /27) ) ; 

g=- ( b / 2 ) + f ; 

h— ■ ( b / 2 ) - f ; 

% Perform necessary manipulations of the coefficients to acquire values 
% in the correct quadrants when taking cube roots of complex quantities 
if (f>0) & ( g<0 ) 

A= ( abs (g) A (1/3) ) *-l; 
else 

A=g A (1/3) ; 

end;%end-if ( f > 0 ) & (g<0) 

if ( f >0 ) & (h<0 ) 

B= (abs (h) A ( 1/3 ) )*-l; 
else 

B=h A (1/3) ; 

end;%end-if (f>0) & (h<0) 

% Calculate and store the three roots of the wave speed ratio squared 
xl=A+B; 

x2= ( - (A+B) /2) + ( ( (A-B) / 2 ) *sqrt (-3) ) ; 
x3= ( - (A+B) / 2 ) — ( ( (A-B) / 2 ) ^sqrt (-3) ) ; 
yl=xl- (p/3 ) ; 
y2=x2- (p/3 ) ; 
y3=x3- (p/3) ; 

speed_ratiosqd_root l=yl ; 
speed_ratiosqd_root2=y2 ; 
speed_rat iosqd_root 3=y3 ; 
speed_rat iosqd_roots= [ y 1 , y2 , y3 ] ; 
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Calculate and store the ratio of longitudinal attenuation coefficient 
squared to the Rayleigh wave number squared 
longsqd_rat iol=l- ( gamma sqd*yl ) ; 
longsqd_ratio2=l- (gammasqd*y2 ) ; 
longsqd_ratio3=l- ( gammasqd*y3 ) ; 

longsqd_rat ios= [ longsqd_ratiol , longsqd_ratio2 , longsqd_ratio3 ] ; 

Calculate and store the ratio of vertical attenuation coefficient 
squared to the Rayleigh wave number squared 
vert sqd_ratio 1=1- yl ; 
vertsqd_ratio2=l-y2 ; 
vertsqd_ratio3=l-y3 ; 

vert sqd_ratios= [ vertsqd_ratiol , vertsqd_ratio2 , vert sqd_ratio3 ] ; 

Perform required logical comparisons to keep only the root that 
satisfies the Rayleigh wave criteria 
suitable= ( longsqd_ratios>=0) & ( vertsqd_ratios>=0 ) ; 
for m=l : 3 

if suitable (m) 

Rayleighsqd_root=speed_ratiosqd_roots (m) ; 
index=m; 

end;%end-if suitable 
end; %end-for m 

Calculate the values of the Rayleigh wave root and the attenuation 
ratios 

Rayleigh_root=sqrt ( Rayleighsqd_root ) ; 

Rayleigh_long_ratio=sqrt ( longsqd_rat ios (index) ) ; 
Rayleigh_vert_ratio=sqrt ( vertsqd_ratios (index) ) ; 

% Calculate the normalized horizontal and vertical particle 
% displacements 

q=Rayleigh_long_ratio; 
s=Rayleigh_vert_ratio; 
l° n g_ con st= (2 + q + s)/((s /v 2)+l); 
vert_const=2/ ( (s A 2 ) +1 ) ; 
depth= [0:0.01:1.6] ; 

surf ace_l a yer= zeros ( 1 , length (depth) ) ; 
krz=depth* (2*pi) ; 

u3= (exp (-q*krz) - ( vert_const *exp (-s*krz ) ) ) ; 

u30= ( l-vert_const ) ; 

u3_u30=u3/u30; 

ul=exp (-q^krz) - ( long_const *exp (-s*krz) ) ; 
ul_u30=ul/ (-q*u30) ; 

% Plot the normalized horizontal and vertical particle displacements 
f igure ( 1 ) ; 

plot (depth, surf ace_layer, depth, u3_u30, depth, ul_u30 ) ; 
xlabel (' depth/Rayleigh wave number'); 
ylabel (’ normalized displacement amplitude'); 
text (0.22,0.05, ' longitudinal ' ) ; text (0.62,0.5, 'vertical ’ ) ; 
title (' Longitudinal and Vertical Particle Displacements'); 
grid on; 
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normalized displacement amplitude 
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APPENDIX A5. TENSILE AND SHEAR STRESSES 



% Filename: rayleighstresses . m 
% Written by: F.E. Gaghan 

% Date Last Modified: 22 December 1997 

% Purpose: This program is used to graphically display the horizontal 
% and vertical particle stresses induced by a propagating 

% Rayleigh wave for a specified value of Poisson's ratio. The 

% curves displayed are the particle stress amplitudes 

normalized by the horizontal particle stress at the surface 
as a function of depth normalized by the Rayleigh wave 
number . 

Initialize required storage vectors 
clear all; 
nu=0 . 25 ; 

gammas qd= (1- (2*nu) ) / (2- (2*nu) ) ; 
g amma s qdi n v= 1 / g amma s qd ; 
lame= (2*nu) / (1- (2*nu) ) ; 

Calculate the coefficients required to solve the Rayleigh wave 
equation which is cubic in the wave speed ratio 

p=-8 ; 

q=24- ( 1 6*gammasqd) ; 
r= ( 16*gammasqd) -16; 
a= ( (3*q)-(p A 2) )/ 3 ; 
b=( (2* (p A 3) )-(9^p*q)+ {27 ^r) ) /27; 
f=sqrt ( ( (b A 2) /4) + ( (a A 3) /27) ) ; 
g=- (b/2 ) +f ; 
h=- (b/2 ) -f ; 

Perform necessary manipulations of the coefficients to acquire values 
in the correct quadrants when taking cube roots of complex quantities 
if (f >0) & (g<0) 

A=(abs(g) A (1/3) ) *-l; 
else 

A=g A (1/3) ; 

end;%end-if ( f > 0 ) & (g<0) 

if ( f >0 ) & ( h<0 ) 

B= ( abs (h) A (1/3) ) *-l; 
else 

B=h A (1/3) ; 

end;%end-if (f>0) & (h<0) 

% Calculate and store the three roots of the wave speed ratio squared 
xl=A+B; 

x2=(-(A+B)/2) + ( ( (A-B) /2) *sqrt (-3) ) ; 
x3= ( - ( A+B ) / 2 ) — ( ( (A-B) / 2 ) *sqrt (-3) ) ; 
yl=xl- (p/3) ; 
y2=x2- (p/3) ; 
y3=x3- (p/3) ; 

speed_ratiosqd_rootl=yl ; 
speed_rat iosqd_root 2=y2 ; 
speed_ratiosqd_root3=y3; 
speed_ratiosqd_roots= [ yl , y2 , y3 ] ; 
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Calculate and store the ratio of longitudinal attenuation coefficient 
squared to the Rayleigh wave number squared 
longsqd_ratiol=l- ( gammasqd*y 1 ) ; 
longsqd_ratio2=l- ( gammas qd*y2 ) ; 
longsqd_ratio3=l- ( gammasqd*y3 ) ; 

longsqd_ratios= [longsqd_ratiol, longsqd_ratio2 , longsqd_ratio3 ] ; 

Calculate and store the ratio of vertical attenuation coefficient 
squared to the Rayleigh wave number squared 
vertsqd_ratiol=l-yl ; 
vertsqd_ratio2=l-y2 ; 
vertsqd_ratio3=l-y3 ; 

vertsqd_ratios= [ vertsqd_ratiol, vertsqd_ratio2 , vertsqd_ratio3 ] ; 

Perform required logical comparisons to keep only the root that 
satisfies the Rayleigh wave criteria 
suitable= ( longsqd_rat ios>=0 ) & ( vertsqd_ratios>=0 ) ; 

for m=l : 3 

if suitable (m) 

Rayleighsqd_root = speed_ratiosqd_roots (m) ; 
index=m; 

end; %end-if suitable 
end; %end-for m 

Calculate the values of the Rayleigh wave root and the attenuation 
ratios 

Rayleigh_root=sqrt ( Rayleighsqd_root ) ; 

Rayleigh_long_rat io=sqrt ( longsqd_rat ios ( index) ) ; 

Rayleigh_vert_rat io^sqrt ( vertsqd_rat ios (index) ) ; 

% Calculate the normalized horizontal and vertical particle stresses 
q=Rayleigh_long_ratio; 
s=Rayleigh_vert_ratio; 
long_const= (2*q*s)/( (s^2) +1) ; 
vert_const=2 / ( (s A 2 ) +1 ) ; 
depth= [0:0.01:1.6] ; 

sur f ace_layer=zeros ( 1 , length (depth) ) ; 
krz-depth* (2*pi) ; 

sigmaxx= ( gamma sqdinv* exp ( -q + krz) ) - 

(lame* (q /N 2) *exp (-q*krz) ) - (2*long_const *exp (-s*krz) ) ; 
sigmaxxO=gammasqdinv- (lame* (q A 2) ) - ( 2* long_const ) ; 
sigmaxx_xxO=sigmaxx/sigmaxxO ; 
sigmazz= ( -gamma sqdinv* (q / "2 ) *exp (-q*krz) ) + 

(lame* exp (-q*krz) ) + ( 2 *long_const *exp (-s*krz) ) ; 
sigma zz_xxO=sigmazz/sigmaxxO ; 

sigmaxz= (-2*q*exp (-q*krz) ) + (2*q*exp (-s*krz) ) ; 
sigmaxz_xxO=sigmaxz/sigmaxxO ; 
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% Plot the normalized horizontal and vertical particle stresses 
figure ( 1 ) ; 

plot (depth, surf ace_layer , depth, sigmaxx_xxO, depth, sigmazz_xxO, 
depth, sigmaxz_xxO) ; 

xlabel { 1 depth/Rayleigh wave number’); 
ylabel ( 'normalized stress amplitude' ) ; 

text (0.06, 0.75, ' longitudinal ’ ) ; text (0.63,0.35, ' shear ' ) ; 
text ( 0 . 42 , 0 . 16, 'vertical'); 

text (0.06,0.7, 'stress' ) ; text (0.63,0.3, ’stress') ; 

text (0 . 42, 0 . 11, ’stress’); 

tit le (' Tensile and Shear Stresses'); 

grid on; 



Tensile and Shear Stresses 
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APPENDIX A6. VARIATION IN RATIO OF LONGITUDINAL TO VERTICAL 
PARTICLE DISPLACEMENT AS A FUNCTION OF POISSON’S RATIO 



Filename: rayleighaxesratio . m 
Written by: F.E. Gaghan 
Date Last Modified: 10 February 1998 
% Purpose: This program is used to graphically display the ratio of the 
% longitudinal to vertical particles displacements at the 

% surface during Rayleigh wave propagation as a function of 

% the Poisson’s ratio. The calculations are for a homogeneous, 

% isotropic, elastic half-space. 

% Establish desired Poisson's ratio values 
clear all; 
nu= [0.0:0.01:0.5] ; 
for i=l : length (nu) 

gammasqd= ( 1- (2*nu (i) ) ) / (2- (2*nu (i) ) ) ; 

% Calculate the coefficients required to solve the Rayleigh wave 
% equation which is cubic in the wave speed ratio 

p=-8; 

q=24 - ( 16*gammasqd) ; 

r = ( 1 6 * gamma s qd ) - 1 6 ; 

a=( ( 3 * q ) - (p~2) ) /3; 

b=( (2* (p^3 ) ) - (9*p*q) + ( 27 *r ) ) /27; 

f =sqrt ( ( (b A 2) /4) + ( (a A 3)/27) ) ; 

g=- (b/2 ) +f ; 

h=- (b/2) -f ; 

% Perform necessary manipulations of the coefficients to acquire 
% values in the correct quadrants when taking cube roots of complex 
% quantities 

if ( f >0 ) & (g<0 ) 

A«(abs(g) A (l/3))*-l; 

else 

A=g~ (1/3)-; 

end;%end-if ( f >0 ) & (g<0) 

if ( f >0) & ( h<0 ) 

B= ( abs (h) " (1/3) ) *-l; 
else 

B=h A (1/3) ; 

end;%end-if ( f >0 ) & (h<0) 

% Calculate and store the three roots of the wave speed ratio 
% squared 
xl=A+B; 

x2= (- (A+B) /2)+ ( ( (A-B) 12) *sqrt (-3) ) ; 
x3= ( - (A+B) / 2 ) — ( ( (A-B) / 2 ) *sqrt (-3) ) ; 
yl=xl- (p/3 ) ; 
y2-x2- (p/3) ; 
y3=x3- (p/3) ; 

speed_rat iosqd_root l=yl ; 
speed_ratiosqd_root2=y2 ; 
speed_rat iosqd_root 3=y3 ; 
speed_rat iosqd_roots= [ yl , y2, y3] ; 
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% Calculate and store the ratio of longitudinal attenuation 
% coefficient squared to the Rayleigh wave number squared 
longsqd_rat iol— 1- ( gammas qd *y 1 ) ; 
longsqd_rat io2=l- (gammasqd*y2 ) ; 
longsqd_rat io3=l- (gammasqd*y3 ) ; 

longsqd_ratios= [ longsqd_rat iol , longsqd_ratio2 , longsqd_rat io3 ] ; 

% Calculate and store the ratio of vertical attenuation coefficient 
% squared to the Rayleigh wave number squared 
vertsqd_rat iol=l-yl ; 
vertsqd_rat io2=l-y2 ; 
vertsqd_ratio3=l-y3 ; 

vertsqd_ratios= [verts qd_r at iol , vertsqd_ratio2 , vertsqd_rat io3 ] ; 

% Perform required logical comparisons to keep only the root that 
% satisfies the Rayleigh wave criteria 

suitable= ( longsqd_ratios>=0 ) & ( vertsqd_ratios>=0 ) ; 

for m=l : 3 

if suitable (m) 

Rayleighsqd_root=speed_ratiosqd_roots (m) ; 
index=m; 

end;%end-if suitable 
end; %end-for m 

% Calculate the values of the Rayleigh wave root and the 
% attenuation ratios 

Rayleigh_root=sqrt ( Rayleighsqd_root ) ; 

Rayleigh_long_ratio=sqrt ( longsqd_ratios (index) ) ; 

Rayleigh_vert_rat io=sqrt ( vert sqd_ratios (index) ) ; 

% Calculate the normalized horizontal and vertical particle 
% displacements and the ratio of horizontal to vertical particle 
% displacements 

q=Rayleigh_long_ratio; 
s=Rayleigh_vert_ratio; 
long_const= (2*q + s)/((s A 2)+l); 
vert_const=2/ ( (s A 2) +1) ; 
u30=l-vert_const ; 
ulO=l-long_const ; 
axes ratio (i)=ul0/(-q*u30) ; 
end; %end-for i 

% Plot the ratio of horizontal to vertical particle displacements as a 
% function of the Poisson's ratio 
figure ( 1 ) ; 
plot (nu, axesratio) ; 
xlabel (' Poisson' s ratio'); 

ylabel (' longitudinal displacement /vertical displacement'); 
tit le (' Longitudinal and Vertical Particle Displacements'); 
axis ( [nu ( 1 ) , nu ( i ) , min ( axesratio) , max (axesratio) ] ) ; 
grid on; 
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Variation in Ratio of Longitudinal to Vertical Particle Displacements 
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APPENDIX A7. VARIATION IN RAYLEIGH WAVES AS A FUNCTION OF 

DEPTH 



% Filename: rayleighwaves3d . m 
% Written by: F.E. Gaghan 
% Date Last Modified: 22 December 1997 

% Purpose: This program produces a 3D graphical representation of 
% Rayleigh waves in an isotropic elastic half space from which 

can be seen the exponential decay with depth into the 
medium . 

% Initialize required space and wave parameters 
clear all; 
xl= [ 0 : 0 . 2 : 8 ] ; 
k= (2*pi*100) / 3 0 0 ; 
kxl=k*xl ; 
w=k*300 ; 
t= [0 : 0 . 2 : 8] ; 
wt=w* t ; 

C=0 . 25 ; 

[XI ] =meshgrid (xl ) ; 

% Loop through various depths in the medium to calculate the horizontal 
% and vertical particle displacements. 

1=0; 

for x3= [ 0 , -0. 125, -0.25, -0. 375] ; 

% Calculate the horizontal particle displacement 
al=exp (0. 8475*k*x3) - (0. 57 7 3* exp (0.3933*k*x3) ) ; 
cl=sin (wt-kxl ) ; 
ul=C* (al *cl ) ; 
ul=ul+x3 ; 

% Calculate the vertical particle displacement 

a3= (-0.8475* exp (0. 8475*k*x3) ) + ( 1 . 4 67 9*exp ( 0 . 3 933* k*x3 ) ) ; 
c3=cos (wt-kxl ) ; 
u3=C* (a3*c3) ; 

% Conditional statement used to display surface layer 
if i==0 

[LX, LY] =meshgrid (xl , u3) ; 

[ L Z ] = 0 * LX ; 
figure ( 1 ) ; 
surf (LX, LY, LZ) ; 
end;%end-if i 

% Plot the waves for each of the four different depths within the 
% medium 

[U3 , U1 ] =meshgrid (ul , u3) ; 

hold on; 

mesh (Xl, U1,U3) ; 

xlabel (' longitudinal axis [m] ' ) ; ylabel (’ transverse axis [m] 1 ) ; 
zlabel (’ depth axis [m] ' ) ; 

title (’ Variation of Rayleigh Waves as a Function of Depth (3D)’); 
view (20,20) ; 
i=i+l ; 

end;%end-for x3 
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depth axis [m] 



Variation of Rayleigh Waves as a Function of Depth (3D) 




transverse axis [m] 



longitudinal axis [m] 



Variation of Rayleigh Waves as a Function of Depth (2D) 
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APPENDIX B. MANUFACTURER SPECIFICATION SHEETS 



The appendix contains the specification sheets from the manufacturer catalogs for 
the various equipment components used during the experimental testing phases. 
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